Equation-of-state dependence of the gravitational-wave signal from the ring-down 

phase of neutron-star mergers 
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Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descrip- 
tions of high-density matter in order to explore the equation-of-state dependence of the postmerger 
ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS 
is the typical outcome of the coalescence of two stars with 1.35 M@ for most candidate EoSs. The 
oscillations of this object imprint a pronounced peak in the gravitational-wave (GW) spectra, which 
is used to characterize the emission for a given model. The peak frequency of this postmerger GW 
signal correlates very well with the radii of nonrotating NSs, and thus allows to constrain the high- 
density EoS by a GW detection. In the case of 1.35-1.35 Mq mergers the peak frequency scales 
particularly well with the radius of a NS with 1.6 Mq, where the maximum deviation from this 
correlation is only 60 meters for fully microphysical EoSs which are compatible with NS observa- 
tions. Combined with the uncertainty in the determination of the peak frequency it appears likely 
that a GW detection can measure the radius of a 1.6 Mq NS with an accuracy of about 100 to 
200 meters. We also uncover relations of the peak frequency with the radii of nonrotating NSs with 
1.35 Mq or 1.8 M@, with the radius or the central energy density of the maximum-mass Tolman- 
Oppenheimer-Volkoff configuration, and with the pressure or sound speed at a fiducial rest-mass 
density of about twice nuclear saturation density. Furthermore, it is found that a determination 
of the dominant postmerger GW frequency can provide an upper limit for the maximum mass of 
nonrotating NSs. The effect of variations of the binary setup are investigated and corresponding 
functional dependences between the peak frequency and radii of nonrotating NSs are derived. With 
higher total binary masses, correlations are tighter for radii of nonrotating NSs with higher masses. 
The prospects for a detection of the postmerger GW signal and a determination of the dominant 
GW frequency are estimated to be in the range of 0.015 to 1.2 events per year with the upcoming 
Advanced LIGO detector. 

PACS numbers: 26.60.Kp, 97.60. Jd, 04.30.Db, 95.85.Sz, 04.25.dk, 95.30.Lz 



I. INTRODUCTION 

Neutron stars (NSs) are the smallest, densest, and 
most compact stellar objects known to exist in the uni- 
verse. The conditions in the cores of NSs at several 
times nuclear saturation density (po = 2.7 ■ 10 14 g/cm 3 ) 
have remained a mystery since the discovery of these ob- 
jects. In particular, the composition and the equation 
of state (EoS) of high-density matter are only incom- 
pletely known [l-3]. This results from the impossibility 
to directly access NS matter by laboratory experiments, 
although some constraints are provided by the properties 
of neutron-rich nuclei and by heavy-ion collisions as 
well as from uncertainties in nuclear forces and in solving 
the strongly interacting many-body problem (see e.g. [l|- 
I)- 

From the astrophysical point of view the EoS is cru- 
cial because it uniquely determines the NS structure and 
therefore the mass-radius relation of NSs, via the stellar 
structure equations (see e.g. HI)- This one-to-one 
correspondence between the EoS and the mass-radius re- 
lation of NSs allows in principle to determine or at least 
to constrain the EoS from the simultaneous measurement 
of the mass and the radius of a NS. While NS masses have 



been measured very precisely in binary systems [J [B-Q] , 
the determination of radii has not yet been achieved with 
high precision (see e.g. [H, [H, [||, Q an d references therein). 
However, mass measurements alone already constrain the 
possible EoSs because all mass-radius relations feature a 
maximum-mass configuration, where more massive ob- 
jects are unstable and collapse to a black hole (BH). The 
detection of a pulsar with (1.97 ± 0.04) M Q via the de- 
termination of the Shapiro delay in a NS-white dwarf 
binary rules out a variety theoretical descriptions of the 

EoS [nj. 

To date about ten NS binaries are known, which con- 
tain at least one star visible as a pulsar. The stellar 
masses in these systems have partially been measured 
very accurately and cluster at about 1.35 M Q [l], 0- 
7] . Because of observational selection effects one expects 
many more NS binaries to exist, which is also predicted 
by theoretical models of the NS population [ll|. Obscr- 
vationally is has been confirmed that the NSs in a binary 
slowly approach each other because of the emission of 
gravitational waves (GWs), which are generated by the 
orbital motion and which carry away angular momentum 
and energy. As a consequence of these losses, it is the fate 
of (isolated) NS binaries to inspiral and finally to merge 
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on timescalcs which arc found to be of the order of several 
hundred Megayears for the known systems [!, H2 • 

NS mergers have been identified as events whose dy- 
namics depend crucially on the high-density EoS (see 
e.g. [13l - [l5j for a review). Consequently, in numerical 
models potentially observable signals of NS coalescences 
have been found to show an EoS dependence. For in- 
stance, short gamma-ray bursts are speculated to origi- 
nate from the BH-torus system forming milliseconds after 
the merging of the binary components [l6j . In particu- 
lar the survival time of the merger remnant before BH 
formation and the torus mass providing the energy in- 
put to the burst has been found to depend on the bi- 
nary masses and the EoS [13, EH > where also the accu- 
rate treatment of thermal effects plays a role [l9|. In 
the most extreme cases, the direct collapse of the merger 
remnant to a BH or a very long phase of stability before 
it becomes gravitationally unstable, might disfavor the 
gamma-ray burst production, a possibility that depends 
sensitively on the mass of the hypermassive remnant rel- 
ative to the NS mass limit of the EoS. It has also been 
shown that the amount of matter becoming gravitation- 
ally unbound from the merger site, for a given binary con- 
figuration is influenced by the adopted model of supranu- 
clear matter [2(j. These ejecta (typically a few 10 -3 M Q ) 
are observationally relevant because nucleosynthetic re- 
actions in the outflow can form the so-called r-process 
elements, heavy neutron-rich nuclei whose astrophysical 
production site has not been identified yet (see [2l[ for 
calculations which take into account in particular the de- 
tailed properties of the EoSs). Observationally not only 
the nucleosynthesis yields are important but also the pos- 
sibility to detect an electromagnetic counterpart of a NS 
merger (2l|-l25j. The signal, which is powered by the 
radioactive decay of the freshly synthesized nuclei, may 
provide information about the host galaxy of an event. 

However, it appears unlikely that the aforementioned 
observations can significantly constrain the properties of 
the high-density matter because the EoS effects are too 
indirect or still not fully explored, and additional infor- 
mation e.g. about the involved binary masses is required. 
More promising in this respect is the detection of GWs 
from inspiraling and merging NS binaries by Advanced 
LIGO (13| or Advanced Virgo [13] going into operation 
within the next years. Indeed, NS mergers are among 
the prime targets for these detectors with expected event 
rates of 0.4 to 400 [28j . It is clear that the determina- 
tion of the involved NS masses will be possible by the 
matched filtering technique [29]. But, it is crucial for 
constraining the EoS to measure simultaneously the mass 
and the radius of a NS. Information on the radius can 
be obtained from the late inspiral signal when finite-size 
effects become important and cause deviations from a 
point-particle behavior and finally terminate the inspiral 
phase (see (30l - [32l | and references therein). In [3l| it was 
worked out that NS radii can be determined with an un- 
certainty of about 1 km. Only symmetric binaries with 
two NSs of 1.35 Mq have been considered there, employ- 



ing simplified EoSs. The effects of the initial mass ratio, 
the total binary mass, and especially of microphysical 
EoSs have been left unexplored. Apart from these open 
questions, alternative methods for determining NS radii 
arc desirable in any case. 

Besides the impact on the inspiral signal, different 
EoSs have been recognized to have a strong influence on 
the postmerger GW signal in the cases when the merg- 
ing stars do not directly form a BH [IH, l33T - l39| . Based on 
calculations for a limited set of EoSs it has been shown 
that the EoS determines the total binary mass for which 
such a prompt collapse does not occur, which may pro- 
vide a handle on the maximum mass of NSs [U, [38j . This 
procedure, however, requires additional knowledge of the 
radius of a 1.4 M Q NS, which is not yet available. Fur- 
thermore, in the simulations of [38| shock heating was 
assumed to be very inefficient. 

In order to identify the EoS influence on the post- 
merger GW emission the full functional dependence of 
the GW signal on the EoS has to be determined. With a 
systematic survey based on a large set of representative 
candidate EoSs we have demonstrated a strong correla- 
tion between the dominant GW frequency of the post- 
merger phase and the radii of nonrotating NSs [4(J . This 
relation allows to deduce the supranuclear EoS when the 
dominant postmerger GW frequency is observationally 
determined, and therefore represents a complementary 
approach to methods relying on the inspiral signal. 

The purpose of this paper is 1) to extend our recent 
study of the EoS effects on the postmerger GW signal by 
considering an enlarged set of 38 microphysical EoSs, 2) 
to investigate the influence of the total binary mass and 
the initial mass ratio, 3) to explore the scatter inherent 
to the relations between the postmerger GW frequency 
and NS radii and thus to specify the accuracy to which 
NS radii can be determined, 4) to evaluate additional 
dependences between the late-time GW signal and EoS 
features, and finally 5) to present details on the employed 
EoSs, the detectability of the GW signature, and the 
involved uncertainties. 

The paper is organized as follows. First we introduce 
the code used for our study focussing in particular on the 
implementation of the EoS. Sect. IIIII outlines details on 
the considered EoSs and the corresponding stellar prop- 
erties of NSs. Then we describe the simulations and the 
analysis of the GW emission including a discussion of the 
detectability (Sect. HV|). In Sect. [V] the relations between 
features of the GW signal and stellar parameters are de- 
scribed. Here, we also discuss the spread in the presented 
relations and interpret the findings. Sect. I VII examines 
direct correlations between the dominant GW frequency 
of the postmerger phase and microphysical properties of 
the EoS. Variations of the binary mass and the mass ratio 
are investigated in Sect. IVlfl Finally, Sect. I Villi presents 
a summary and our conclusions. 

Note that in this paper we refer to the gravitational 
mass of an object in isolation if not specified otherwise. 
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II. METHODS AND INITIAL DATA 

The simulations of our survey are conducted with a 
three-dimensional relativistic Smoothed Particle Hydro- 
dynamics (SPH) code, where the gravitational fields are 
obtained by solving the Einstein eq uations within the 
conformal flatness approximation [4lll43l . A detailed de- 
scription of the code can be found in [2(J [53, E3] . Within 
the SPH approach the stellar fluid is modelled by a set 
of particles with constant rest mass, which are consid- 
ered not to be pointlike but to be spread out over a 
small spatial domain. Hydrodynamical quantities are as- 
signed to the particles and are evolved comoving with 
the fluid, i.e. according to the Lagrangian formulation of 
hydrodynamics. Specifically, we evolve a set of so-called 
conserved quantities, the conserved rest- mass density p*, 
the conserved specific momenta Uj and the conserved en- 
ergy density r, by a system of coupled ordinary differen- 
tial equations. These conserved quantities are functions 
of the metric potentials and the "primitive" quantities, 
the rest-mass density p, the coordinate velocities Vi, the 
specific internal energy e, and the pressure P. In ev- 
ery timestep the primitive variables are recovered from 
the evolved quantities (p*,Ui,r) by an inversion proce- 
dure, which requires to invoke the EoS P = P(p,Y c ,e) 
with the electron fraction Y e to close the system of hy- 
drodynamical equations. This step is necessary because 
the pressure explicitly occurs in the evolution equations. 
The electron fraction is evolved according to = 0, i.e. 
advected with the fluid, because the timcscale of lepton 
number transport through neutrinos is long compared to 
the dynamical timescale (at least for the bulk of the stel- 
lar matter). 

Concerning technical aspects the EoSs employed in our 
survey can be grouped into four categories. 

(i) Some EoSs provide the full dependence P = 
P{p,Y c ,T), where the temperature T is obtained 
by inverting e = e(p, Y e ,T). These EoSs are avail- 
able in the form of tables, where we use a trilinear 
interpolation scheme. 

(ii) Many published EoSs describe high-density mat- 
ter at zero temperature and in equilibrium with 
respect to weak interactions (so-called neutrinolcss 
beta-equilibrium with the neutrino chemical poten- 
tial jU w = 0). Thus, for these "barotropic" EoSs the 
pressure is given as a unique function of the rest- 
mass density or, alternatively, the energy density 
e = p(l + e). In these cases a thermal component of 
the EoS can be taken into account approximativcly 
by means of an ideal gas-like ansatz. This is nec- 
essary because temperature effects become impor- 
tant when the NSs collide. To this end we employ a 
scheme detailed in [f9j], where a thermal contribu- 
tion of the specific internal energy is computed via 
£th = e — £coid(p)- The rest-mass density p and the 
specific internal energy e are given by the evolution 
of the hydrodynamical scheme. The cold contri- 



bution e co id is determined by p only by calling the 
"cold" barotropic EoS. The thermal contribution 
to the pressure reads P t h = (Pth — 1 )p^th with the 
ideal-gas index r t h- It is added to the cold compo- 
nent to yield the total pressure P = P co \d{p) + Pth- 
Here, the cold part for a given rest-mass density is 
again obtained only from the barotropic EoS. The 
validity and the involved uncertainties of this ap- 
proach have been discussed in There it has 
been shown that in the context of NS mergers a 
value of T t h = 2 represents a suitable choice to 
mimic thermal effects of nuclear matter. We there- 
fore use Tth = 2 if not mentioned otherwise. In par- 
ticular, the collapse timescale of the merger rem- 
nant is reproduced very well by this choice of Tth- 

(iii) Additional barotropic EoSs which are not avail- 
able to us as tables can be included by their 
parametrized versions provided in p3 | in form of 
piecewise polytropes. Also for this class of mod- 
els we make use of the thermal ideal-gas extension 
outlined above. One should bear in mind that the 
polytrope fits imply approximations that will have 
to be considered later in our discussion. First, the 
parametrized versions do not perfectly match the 
underlying microphysical EoSs (see [44| ) and hence 
simulation results may differ because of the use of 
the piecewise polytrope instead of the original EoS. 
Second, the properties (e.g. radii) of cold NSs ob- 
tained from the fitted versions are slightly different, 
thus shifting our results when presented in depen- 
dence of these NS properties. 

(iv) The fourth class of EoSs is based on the calcula- 
tions of [If| . Detailed computations of the EoS be- 
low nuclear saturation density using chiral effective 
field theory interactions are extended by piecewise 
polytropes for higher densities requiring a maxi- 
mum mass of NSs above 1.97 M Q . At densities 
below 5 • I0 10 g/cm 3 we extend the original tables 
with an additional crust. The resulting barotropic 
relations are again supplemented with the ideal-gas 
component for thermal effects. 

Note that the classes (ii) to (iv) involve an additional 
approximation beside the treatment of thermal effects 
because matter is assumed to adjust instantaneously to 
beta-equilibrium as given by the barotropic EoSs (see the 
discussion in [H|). In these cases the electron fraction is 
not used as an input for the EoS call because the pressure 
is entirely determined by the rest-mass density and the 
specific internal energy. 

The details of the different EoSs employed here, their 
properties and the underlying microphysical models are 
discussed in Sect. IIIII Note that for indicating the dif- 
ferent classes of EoSs in figures we use the same color 
scheme throughout this paper (class (i) in red, class (ii) 
in black, class (iii) in green, class (iv) in blue). 

General relativistic gravity cannot be treated within 
the SPH scheme. Therefore, the five coupled nonlinear 
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elliptic equations resulting from the Einstein equations 
by imposing a confromally flat spatial metric are solved 
on an overlaid grid covering well the binary orbit and 
the postmerger remnant at later times. By the confor- 
mal flatness approximation the effects of GWs are ne- 
glected, for which reason an additional GW backreaction 
scheme is implemented in the code to account for the 
loss of energy and angular momentum by the emission of 
gravitational radiation. This procedure requires to solve 
additional elliptic equations to compute corrections to 
the conformally flat metric. Details are reported in (2pj . 

Magnetic fields or neutrino transport physics are not 
included in our model, because it has been shown that 
these effects have only little influence on the dynamics 
of a merger as long as the initial field strength inside 
the stars is below fO 14 G and only a fe w dy namical 
timescales are tracked after the merging [4a - l50 |. 

Our simulations start from a quasi-equilibrium orbit 
about two to three revolutions before the actual merging 
of the binary components. Initially we impose a nonro- 
tating (irrotational) velocity profile of the NSs, because 
the inspiral time is not sufficient to enforce synchronous 
rotation (5f], [HJ ■ This state can be considered to be a 
good approximation also for NSs spinning with millisec- 
ond periods because even in such cases the rotation is still 
slow compared to the orbital period. Furthermore, the 
initial stars are assumed to have zero temperature and 
to be in neutrinoless beta-equilibrium. The two NSs are 
set up typically with about 340,000 SPH particles, but 
we perform additional higher resolved calculations to en- 
sure that our results do not depend on the resolution (see 
Subsect.IVB]). 



III. EQUATIONS OF STATE AND STELLAR 
PROPERTIES 

The EoSs used in our survey are listed in Tab.|H where 
it is also indicated to which technical category a partic- 
ular EoS belongs. Seven out of the 38 models include 
thermal effects consistently (class i). The remaining EoSs 
have been computed for high-density matter at zero tem- 
perature and in equilibrium with respect to weak inter- 
actions (barotropic EoSs). 

The EoSs of classes (i), (ii) and (iii) have been de- 
rived within different frameworks, including nonrelativis- 
tic and relativistic nuclear energy-density functionals, 
Brueckncr-Hartrcc-Fock calculations, and phenomeno- 
logical models, such as the liquid-drop model (see e.g. 0, 
Sll for an overview) . All except for two microphysical 
models describe interacting neutrons, protons and elec- 
trons (possibly also muons) , where some also include hy- 
perons (Glendnh3, H3, H4). MIT60 and MIT40 consider 
strange quark matter and are discussed below. In addi- 
tion, the nonzero temperature EoSs include positrons and 
photons. Also the models for the (nuclear) interactions 
have been varied, which implies that the EoSs explore a 
range of nuclear experimental data. Some EoSs take into 



account three-body forces. 

Figure [T] shows the 6 EoSs of class (iv). The dark 
blue band at low energy densities e shows the uncer- 
tainty band of the pressure obtained from microscopic 
calculations based on chiral two- and three-nucleon in- 
teractions [i^, [54| . This band is extended in a general 
way to higher densities using piecewise polytropes as a 
function of the rest- mass density [45|, [55[. The exten- 
sions are constrained by the causality of the polytropes 
and the heaviest observed NS [l(|. The light blue ex- 
tension band covers all possible EoSs that are compat- 
ible with the constraints based on chiral effective field 
theory interactions at low densities and are able to sup- 
port NSs of mass M ^ 1.97M© [H, HU- In addition, 
we overlay a dark blue extension band, which is obtained 
by requiring a maximum NS mass M 2.4Af Q . The 6 
EoSs of class (iv) within these bands are a representa- 
tive sample of EoSs characterized by different stiffness 
fulfilling these constraints. For our choice of polytopcs 
the EoSs become acausal (the speed of sound exceeds 
the velocity of light) beyond a density pumit where the 
speed of sound is w s (/Cii m it) = c. We replace this unphys- 
ical regime by P = e — const., which is by construction 
causal, and choose the constant in a way that the pres- 
sure is continuous. The revised parts of the EoSs can be 
recognized by the short pieces of the lines sticking out 
of the blue bands at high densities. (The upper edge 
of the blue bands (parallel to the dotted P = e line) 
indicates pu m it-) Note that the EoSs are only replaced 
at very high densities slightly below the central density 
Pmax of the maximum-mass configuration of static NSs 
shown by the black dots. Because of this, the density 
regime above piimit does not play a role during the evo- 
lution of NS mergers for the cases considered here, and 
the replacements of the acausal part of the EoSs affect 
only the exact values of stellar properties of very massive 
static NSs. 

The thermodynamic properties of all EoSs available for 
our study are illustrated in Fig. [5] showing the pressure 
P as a function of the rest-mass density p. One observes 
that in the range of one to several times nuclear satu- 
ration density (p = 2.7 • I0 14 g/cm 3 ) the spread of the 
pressure between the various models is nearly one order 
of magnitude. It has been pointed out that the pressure 
at a density between one and two times po serves as a 
characteristic quantity that correlates with the radius of 
a 1.4 M e NS JJ[H3|. We provide the pressure at 1.85p in 
Tab. HI which is chosen in order to be consistent with the 
fiducial rest-mass density for the piecewise polytrope fits 
presented in (44|. It ranges from 1.18 TO 34 dyne/cm to 
9.56 -10 34 dyne/cm . In comparison to the high-density 
regime the pressure differences decrease towards lower 
densities (not shown in Fig. [2]), see e.g. pi. l45j|. 

As another EoS property we display the speed of sound 
in Fig. [3] for all EoSs. The diversity among the models is 
obvious, where the differences are larger at higher densi- 
ties. We extract the sound speed at 1.85po a s a charac- 
teristic quantity for all EoSs and give the values in Tab.fl] 
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FIG. 1: Pressure versus energy density for the 6 EoSs of class 
(iv) (red lines). The dark blue band at low densities shows 
the results of microscopic calculations based on chiral effective 
field theory interactions [45|, [54|. The light blue extension 
band covers possible polytropic extensions compatible with 
causality and NSs of M ^ 1.97M0. The darker blue extension 
is for M ^ 2.4M Q [H, El- The black dots mark the central 
density of the maximum-mass configuration and the dotted 
line indicates the causal limit P = e. 
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FIG. 2: Pressure as a function of the rest-mass density for 
all candidate EoSs at zero temperature and in neutrinoless 
beta-equilibrium. Microphysical EoSs which provide the full 
temperature dependence are shown in red (class i), microphys- 
ical EoSs that are given at zero temperature and in neutri- 
noless beta-equilibrium are displayed in black (class ii), green 
lines correspond to piecewise polytropes (class iii) and EoSs 
of class iv are depicted in blue. Dashed lines are used for 
strange quark matter EoSs. The vertical line marks a rest- 
mass density of p = 10 14 ' 7 g/cm 3 = 1.85po- The filled circles 
indicate the central rest-mass density of the maximum-mass 
TOV solution for every EoS. 
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FIG. 3: Sound speed for all EoSs used in our survey. Colors 
and line style have the same meaning as in Fig. [2] The vertical 
line marks a rest-mass density of p = 10 14 ' 7 g/cm 3 = 1.85po- 
The filled circles indicate the central rest-mass density of the 
maximum- mass TOV solution for every EoS, if not visible the 
EoS becomes acausal for stable TOV solutions. 



The values vary between 0.29 and 0.70 times the speed of 
light. Note that for some EoS models the sound velocity 
exceeds the speed of light at high densities, i.e. the EoSs 
become acausal (e.g. the eosAU, eosUU and APR EoS). 

The variety of the microphysical models employed 
in our study is also reflected in the mass-radius rela- 
tions of non-rotating NSs, which are obtained by solv- 
ing the relativistic equations of hydrostatic equilibrium 
(Tolman-Oppenhcimcr-Volkoff (TOV) equations). The 
TOV equations yield a unique relation between the NS 
mass and radius for each EoS, which allows to character- 
ize an EoS by the resulting stellar properties. Figure [4] 
displays the TOV solutions for all EoSs. Red curves cor- 
respond to EoSs whose temperature dependence is given 
by the microphysical model, while fully microphysical 
zero-temperature EoSs (which in the simulations are sup- 
plemented with an ideal-gas component to mimic thermal 
effects) are shown with black lines. The green curves ex- 
hibit the mass-radius relations of EoSs of class (iii), while 
EoSs belonging to class (iv) are shown in blue. 

All mass-radius relations feature a maximum-mass 
configuration. More massive objects collapse to a BH 
if not additionally supported, e.g. by rotation or ther- 
mal effects. The corresponding solutions are marked with 
symbols in Fig. 21 and their stellar properties are listed in 
Tab. [TJ Furthermore, Tab. U gives the central energy den- 
sity of these maximum-mass configurations, they range 
from 1.35 • 10 15 g/cm 3 to 3.80 • 10 15 g/cm 3 . The cor- 
responding rest-mass densities are indicated by circles in 
Figs. [5] and [H Note that EoSs where the dot is not visible 
in Fig. [3] become acausal (speed of sound larger than the 
speed of light) at densities below the maximum density. 

Our survey considers EoSs that lead to maximum 
masses in the range of 1.79 Mq to 3.00 M Q and to corre- 
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FIG. 4: Mass-radius relations for all EoSs with the gravita- 
tional mass M in isolation and the areal radius R. The color 
scheme is the same as in Fig. [2] The dashed lines denote 
mass-radius relations for strange quark matter EoSs. The 
horizontal line corresponds to the observed 1.97 M@ NS [Tol | . 
For EoSs where the merger of two stars with 1.35 Mq leads 
to the prompt formation of a BH the maximum-mass configu- 
ration is indicated by a cross. Maximum-mass configurations 
depicted by a circle correspond to EoSs where in the simu- 
lation of this binary setup the formation of a differentially 
rotating object is found. 



sponding radii, denoted as i? max , spanning from 8.65 km 
to 14.30 km (Fig. 0] Tab. There has not been any 
special selection procedure for the EoSs, except that we 
require a maximum mass larger than roughly 1.8 Mq. 
This being fulfilled we include every EoS that is available 
to us. The lower bound of about 1.8 Mq is motivated 
by the discovery of a NS with a gravitational mass of 
(1.97 ±0.04) M Q [k|. This measured mass is indicated 
as horizontal line in Fig. This detection practically 
rules out some EoSs of our sample with M max below the 
limit. We do not dismiss such excluded models because 
they may still provide a viable model at lower densities 
(see also Sect. IV Cp . For instance, during the first 5 ms 
after merging the central density in the merger remnant 
described by the excluded LS180 EoS remains below the 
central density of a nonrotating 1.5 Mq NS modeled by 
this EoS. For such "low- mass" stars the mass-radius re- 
lations of excluded EoSs arc partially similar to those ob- 
tained from EoSs compatible with the observation of [l(| ■ 
Hence, in the corresponding density regimes relevant for 
the low-mass stars and the merger remnant such EoSs 
can still yield a viable description of high-density matter. 
In addition to that, the inclusion of EoSs with relatively 
low M max extends (maybe artificially) the range of varia- 
tions of stellar parameters, and correlations between NS 
properties and GW characteristics that hold over a wider 
parameter range can be inferred easier. We note that all 



of the four technical EoS categories cover a similar range 
of stellar parameter values. Only the mass-radius rela- 
tions of class (iv) lie in a more narrow band, which was 
the main result of [I|| . 

A common feature of most EoSs is a relatively small 
variation of the NS radius between about 0.5 Mq and 
about (M max — 0.5 Mq). This suggests to use the radii 
in this mass range as a characteristic feature of a given 
EoS. 

Finally, the MIT60 and MIT40 EoSs deserve a com- 
ment. These models describe absolutely stable strange 
quark matter within the MIT bag model [13, HH, i.e. a 
deconfined quark phase with an energy per baryon lower 
than the one of nucleonic matter (E/A =860 MeV for 
MIT60 and E/A =844 MeV for MIT40). As a conse- 
quence of the strange matter hypothesis [59l IfjOl ] under- 
lying these two EoSs, the compact stars observed in the 
universe, commonly referred to as NSs, would actually be 
strange quark stars (consisting of strange quark matter) . 
This possibility has not yet been ruled out theoretically 
or observationally (see e.g. [H, H, [53[ for details and for 
observational consequences discriminating this scenario 
from ordinary NS; see [13, HH for the consequences of this 
hypothetical state of matter in the context of compact bi- 
nary mergers). As a striking difference to nucleonic NSs, 
strange quark stars show an inverse mass-radius relation 
typical of this class of objects because of the self-binding 
of strange quark matter. The particular model MIT60 
with M max = 1.88 Mq is excluded by the observation of 
the two-solar-mass pulsar. The MIT40 EoS, however, is 
compatible with present knowledge. For the MIT40 EoS 
belonging to class (ii), we adopt r t h = 1.34. 

Note that throughout this paper we use the more com- 
mon term NS instead of compact star for all compact stel- 
lar objects including strange quark stars. With "purely" 
or "fully" microphysical EoSs we refer to models of class 
(i) or (ii), which do not involve piecewise polytropes (see 
Sect. HI]) . Moreover, in this paper "accepted" EoSs de- 
note models which are compatible with the detection of 
the 1.97 Mq NSs taking into account the error bars of 
the observation by fioj . 



IV. SIMULATIONS 

A. Dynamics 

According to pulsar observations 0, 0] and population 
synthesis studies [ll| binaries of two NSs with a gravita- 
tional mass of about 1.35 Mq each are the most abun- 
dant systems in the binary NS population. Therefore, we 
choose a symmetric binary with Mi = M^ = 1.35 Mq 
and simulate for all EoSs discussed in Sect. IIIII the late 
inspiral phase, the merging, and the early postmerger 
evolution of this system until an approximately station- 
ary state has formed (10 to 20 ms after merging). The 
inspiral is driven by the loss of angular momentum and 
energy due to the GW emission and lasts between some 
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TABLE I: Models for the EoS with their thermodynamical properties and resulting stellar parameters. M max is the maximum 
mass of nonrotating NSs with the radius i? max corresponding to this maximum-mass configuration. The next three columns 
provide the radii of nonrotating NSs with 1.35 Mq, 1.6 Mq and 1.8 Mq. P is the pressure at a fiducial rest-mass density of 
1.85 times nuclear saturation density. The sound speed at the same reference density is given in the eighth column. The central 
energy density of the maximum- mass configuration of the TOV solutions is listed as e max . "Class" identifies how an EoS is 
implemented in our code (see Sect. 11111 for details). References for the EoSs are listed in the last column. For EoSs marked 
with an "x" in parentheses the merging of two NSs with 1.35 Mq leads to the direct formation of a BH. 



EoS 


M max 
[M @ ] 


-^max 

[km] 


-Rl.35 

[km] 


-Rl.6 

[km] 


-Rl.8 

[km] 


P(1.85p ) 

r-i /-\.'-S4 l / 2i 

[10 dyne/ cm ] 


« s (1.85po) 
[c] 


^max 

[10 g/cm J ] 


class 


Ref. 


GS1 


2.75 


13.27 


14.72 


14.79 


14.81 


6.9920 


0.6348 


1.55 


i 


[62] 


LS375 


2.71 


12.34 


13.56 


13.71 


13.78 


5.3904 


0.5791 


1.75 


i 


[63] 


Shen 


2.24 


12.63 


14.64 


14.53 


14.35 


5.3480 


0.5202 


1.40 


i 


[64] 


GS2 


2.09 


11.78 


13.38 


13.31 


13.13 


4.2195 


0.4847 


2.04 


i 


[65] 


LS220 


2.04 


10.61 


12.64 


12.43 


12.10 


3.3569 


0.4170 


2.54 


i 


[63] 


MIT60 


1.88 


10.35 


10.43 


10.74 


10.75 


1.4654 


0.5630 


2.20 


i 


[66] 


LS180 


1.83 


10.04 


12.13 


11.65 


10.68 


2.9883 


0.3840 


2.70 


i 


[63] 


eosL 


2.76 


14.30 


15.74 


15.77 


15.79 


9.5635 


0.5900 


1.35 


ii 


[67] 


eosO 


2.39 


11.56 


12.85 


12.87 


12.83 


3.9106 


0.4786 


2.04 


ii 


[68] 


MIT40 


2.28 


12.56 


12.08 


12.59 


12.89 


5.0864 


0.5687 


1.49 


ii 


[69, 70] 


cosUU 


2.21 


9.84 


11.18 


11.12 


11.00 


1.8118 


0.3477 


2.80 


ii 


[71] 


SKA 


2.21 


11.17 


13.07 


12.91 


12.68 


3.5929 


0.4446 


2.30 


ii a 


[63] 


APR 


2.19 


9.90 


11.33 


11.25 


11.12 


1.9152 


0.3533 


2.80 


ii 


[72] 


eosAU (x) 


2.14 


9.45 


10.44 


10.42 


10.33 


1.1812 


0.2906 


3.00 


ii 


[71] 


BurgioNN 


2.07 


10.58 


11.99 


11.87 


11.67 


2.4720 


0.3936 


2.50 


ii 


[73] 


Sly4 


2.05 


10.01 


11.79 


11.59 


11.29 


2.4012 


0.3830 


2.85 


ii 


[74] 


Glendnh3 


1.96 


11.48 


14.52 


13.98 


13.22 


4.4281 


0.4326 


2.35 


ii 


[75] 


BBB2 


1.92 


9.55 


11.30 


10.97 


10.49 


2.0949 


0.3549 


3.15 


ii 


[76] 


eosC 


1.87 


9.89 


12.06 


11.57 


10.79 


2.6024 


0.3693 


3.04 


ii 


[77] 


eosWS (x) 


1.85 


9.58 


10.99 


10.72 


10.16 


1.9070 


0.3357 


3.09 


ii 


[71] 


FPS (x) 


1.80 


9.30 


10.91 


10.52 


9.38 


1.8720 


0.3344 


3.35 


ii 


[78] 


MSI 


2.77 


13.39 


14.99 


15.04 


15.04 


7.1495 


0.6218 


1.60 


iii 


[44] 


MSlb 


2.76 


13.30 


14.59 


14.69 


14.74 


7.0970 


0.6194 


1.60 


iii 


[44] 


MPA1 


2.47 


11.35 


12.49 


12.54 


12.54 


3.0980 


0.4693 


2.10 


iii 


[44] 


APR3 


2.38 


10.73 


12.04 


12.04 


12.00 


2.4451 


0.4202 


2.35 


iii 


[44] 


ENG 


2.25 


10.41 


12.05 


11.99 


11.87 


2.7104 


0.4136 


2.55 


iii 


[44] 


H4 


2.02 


11.72 


13.95 


13.74 


13.35 


4.6293 


0.4447 


2.15 


iii 


[44] 


ALF2 


1.98 


11.48 


12.78 


12.72 


12.49 


4.0893 


0.4396 


2.10 


iii 


[441 


ALF4 


1.94 


10.79 


11.60 


11.55 


11.43 


2.0437 


0.3789 


2.30 


iii 


[44] 


MS2 


1.80 


11.81 


14.25 


13.52 


12.10 


3.9987 


0.4088 


2.20 


iii 


[44] 


APR2 (x) 


1.79 


8.65 


10.13 


9.74 




1.3266 


0.2887 


3.80 


iii 


[44] 


H3 


1.79 


11.77 


13.95 


13.48 




4.3852 


0.4153 


2.15 


iii 


[44] 


Heb6 


3.00 


13.17 


13.33 


13.54 


13.69 


5.7020 


0.6967 


1.50 


iv 


[45] 


Heb5 


2.79 


12.25 


12.38 


12.59 


12.73 


3.3927 


0.5537 


1.70 


iv 


[45] 


Heb4 


2.74 


12.08 


12.51 


12.65 


12.74 


3.0665 


0.5270 


1.80 


iv 


[45] 


Heb3 


2.42 


10.94 


12.03 


12.09 


12.09 


2.7306 


0.4414 


2.25 


iv 


[45] 


Heb2 


2.22 


10.15 


11.42 


11.35 


11.24 


1.8099 


0.3625 


2.65 


iv 


[45] 


Hebl 


2.13 


9.71 


10.81 


10.77 


10.66 


1.4685 


0.3289 


2.90 


iv 


[45] 



a Thc full temperature dependence of this EoS is available for T < 
30 MeV, which is exceeded in NS merger simulations. Therefore, 
we use the ideal-gas ansatz for thermal effects combined with the 
zero-temperature slice of this EoS. 
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ten and several hundred million years depending on the 
initial binary separation. Only shortly before the merg- 
ing takes place, when the orbital period has reached mil- 
liseconds, finite-size and post-Newtonian effects become 
important. 

It is known that the remnant formation after the co- 
alescence depends on the total binary mass M to t = 
Ml + M 2 and the EoS (see e.g. [H[). In this and the fol- 
lowing two sections the total binary mass is fixed, thus 
the outcome depends on the EoS only. We distinguish 
two scenarios. For four out of the 38 EoSs considered 
in our study, we observe the direct formation of a BH 
within less than a ms after the stars come in contact, be- 
cause the merged object cannot be stabilized against the 
gravitational collapse. This "prompt collapse" scenario 
is found only for EoSs that lead to relatively small R max 
(EoSs labelled with an "x" in Tab. UJ). The maximum- 
mass configurations of these EoSs are marked with a cross 
in Fig. |U 

In the simulations with the remaining EoSs the col- 
liding stars form a differentially rotating object (DRO), 
which is supported against the gravitational collapse - at 
least for a longer period - by centrifugal forces. Further- 
more, shock heating at the contact interface and com- 
pression lead to an increase of the temperature to sev- 
eral 10 MeV (in some cases even more than 100 McV), 
which has an additional stabilizing effect. After angular 
momentum is partially lost by GW emission and redis- 
tributed to the outer remnant layers, the merger rem- 
nant may collapse to a BH. This "delayed collapse" hap- 
pens typically after several 10-100 ms and only takes 
place if the mass of the DRO exceeds the upper mass 
limit of rigidly rotating NSs, which for most EoSs is 
about 1.2 M max @, H, [73. We explicitly point out that 
for some EoSs the total remnant mass may not overstep 
this threshold and the DRO finally will settle to a uni- 
formly and rapidly rotating NS. In these cases the 1.35- 
1.35 Mq NS binaries might be excluded as progenitors of 
short gamma-ray bursts, at least for models of gamma- 
ray bursts that rely on the formation of a BH with a hot 
accretion torus (see [IB] for a review). Rigidly spinning 
"supermassive" NSs (whose gravitational mass is larger 
than M max ) may form BHs on much longer timescales 
due to the loss of angular momentum e.g. by electromag- 
netic emission 2, 5] . Note that the gravitational mass of 
a remnant resulting from a 1.35-1.35 Mq binary is of the 
order of 2.6 Mq because the gravitational mass of the 
initial binary is already below 2.7 Mq. Moreover, energy 
is radiated away by GWs, matter becomes gravitation- 
ally unbound during merging, and energy converted into 
heat is lost by neutrino emission. The latter effect is not 
modelled here. 

Independent of whether a BH is the final outcome, 
NS mergers evolve dynamically similarly for all systems 
which form a DRO. As the stars approach each other 
increasingly faster during the late inspiral phase, the 
binary components get more strongly deformed and fi- 
nally collide with a relatively big impact parameter. The 




5 10 15 20 25 



t [ms] 

FIG. 5: Gravitational-wave amplitude of the plus polarization 
measured along the polar axis at a distance of 20 Mpc for the 
simulation with the Shen EoS. 



bulk matter of the initial stars assembles into a rotat- 
ing double-core structure, where the dense cores bounce 
against each other. While matter is shed off from the 
surface to feed a halo around the central object, the 
two cores merge into a single core after a few bounces. 
On a timescale of several milliseconds the oscillations of 
the initially highly deformed remnant are damped and 
an approximately stationary, axisymmetric object is left, 
which is still rotating differentially and ringing with much 
lower amplitude. We simulate all models resulting in a 
DRO for about 15 ms after the plunge until the nearly 
stationary phase is reached. In none of these calculations 
we actually observe the delayed collapse except for the 
MIT60 EoS, where the merger remnant forms a BH after 
about 3 ms, and for the H3 EoS, where the DRO is stable 
for 25 ms. 



B. Analysis of the gravitational- wave signal 

The dynamical stages of a merger can also be identi- 
fied in the GW signal, which is computed by means of 
a modified quadrupole formula that takes into account 
post-Newtonian effects [20| • A typical waveform is shown 
in Fig. [5] for the Shen EoS as measured perpendicular to 
the orbital plane at a distance of 20 Mpc. Fig. [5] displays 
the plus polarization, the cross polarization looks very 
similar but is phase-shifted by 7r/2. The characteristic in- 
spiral phase with an increasing frequency and an increas- 
ing amplitude until 5 ms is clearly visible. Following the 
plunge the ringdown of the postmcrger remnant can be 
seen by high-frequency oscillations that cease over f to 
15 ms until the object has reached approximate stationar- 
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FIG. 6: Gravitational-wave amplitude of the plus polarization 
measured along the polar axis at a distance of 20 Mpc for the 
simulation with the Sly4 EoS. 
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f [kHz] 



FIG. 7: Fourier spectrum of the plus polarization of the GW 
signal at a distance of 20 Mpc for the Shen EoS. The thick 
line displays the spectrum computed from the signal of the 
full simulation time, while the thin line shows the spectrum 
of the postmerger phase only. The dashed lines give the unity 
SNR sensitivity curve of Advanced LIGO (red) and of the 
Einstein Telescope (black). 



ity. For rather stiff EoSs with R max > 11 km we observe 
a low-frequency modulation of the postmerger signal as 
in Fig. [5] In contrast, for soft EoSs with relatively small 
-Rmax such a feature is absent or less pronounced as for 
instance in the waveform calculated for the model with 
the Sly4 EoS, which is shown in Fig. [5] The reason for 
this difference could not be clarified and deserves further 
investigation, e.g. in terms of an oscillation mode analy- 
sis as in [8(| • Furthermore, one recognizes in Fig. [S] that 
for this EoS the inspiral phase lasts about 2 ms longer 




f [kHz] 



FIG. 8: Fourier spectrum of the plus polarization of the GW 
signal at a distance of 20 Mpc for the Sly4 EoS. The thick 
line displays the spectrum computed from the signal of the 
full simulation time, while the thin line shows the spectrum 
of the postmerger phase only. The spectrum of a simulation 
starting 5.5 revolutions before merging is given by the blue 
line. The dashed lines give the unity SNR sensitivity curve of 
Advanced LIGO (red) and of the Einstein Telescope (black). 

and ends with a larger wave amplitude, which is a conse- 
quence of the higher compactness of the inspiraling NSs. 
Note that all simulations start with the same initial co- 
ordinate distance between the stars. 

To obtain spectral information of the GW signal we 
compute the dimensionlcss quantity /i ff.+ = h + (f)f 
with the Fourier transformed waveform h + (f) of the plus 
polarization. The results are displayed in Fig. [7] and 
Fig. [5] for the Shen EoS and the Sly4 EoS, respectively. 
The thick lines show the spectra calculated from the 
whole signal during the simulation time, while the thin 
lines correspond to the postmerger phase alone. (Note 
that slight differences between the spectrum for the post- 
merger signal and the spectrum for the full signal at the 
highest frequencies depend on the time chosen to define 
the beginning of the postmerger phase.) The dashed lines 
belong to the unity signal-to-noise ratio (SNR) sensitivity 
curves for Advanced LIGO (red) (broadband configura- 
tion) [2(| and the planned Einstein Telescope (black) [U . 
Note that the low-frequency part of the spectra below 
1 kHz computed for the full signal is not reliable be- 
cause our simulations start only a few revolutions before 
merging and therefore lower frequencies of the preceding 
inspiral phase arc underreprescntcd. 

A pronounced peak in the spectra at a frequency / pea k 
between about 2 and 4 kHz is generic to all models re- 
sulting in a DRO. It has a frequency of / pca k = 2.19 kHz 
for the Shen model and of / pca k = 3.32 kHz for the Sly4 
EoS. This feature clearly originates from the postmerger 
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phase and is produced by the violent oscillations of the 
deformed merger remnant. As can be seen in the spectra 
of the postmerger evolution, the structure at / poa k is ac- 
companied by additional, weaker peaks. In a first step to- 
wards astroseismology of NS merger remnants it has been 
shown that the remnant can be considered as an isolated 
selfgravitating object, where certain oscillation modes 
are excited by the formation process of the object [80j |. 
The different structures in the postmerger spectrum are 
caused by the various oscillation modes in the merger 
remnant and their couplings. The study of has also 
identified the mode corresponding to the peak frequency 
/peak as the fundamental quadrupolar fluid mode. 

In Tab.lTTlthe peak frequencies for all models forming a 
DRO are given, where the lowest peak frequency is found 
for the eosL model with / poa k = 1.84 kHz. The highest 
peak frequency of 3.73 kHz is obtained in the simulation 
with the BBB2 EoS. In addition, the full width at half 
maximum (FWHM) is provided for every model. Typi- 
cally, the FWHM is in the range of 100 to 200 Hz, with 
the exception of the MIT60 calculation yielding a broad 
peak with a FWHM of 460 Hz. 

For the Shen EoS we have performed higher resolved 
simulations with 550,000 and 1,270,000 SPH particles, 
which confirm that / pca k is determined to an accuracy of 
about one per cent (see also the discussion and figure 2 
in [4(| and [HI). Moreover, a run starting 3.5 orbits be- 
fore the merging of the binary components does not yield 
differences in the peak frequency (see [Hj ) ■ Additionally 
we show in Fig. [8] the resulting spectrum of a calculation 
starting 5.5 orbits before merging for the Sly4 EoS (blue 
line). The pcak frequency is unaffected by the duration 
of the simulated inspiral phase (difference in / pca k below 
0.5 per cent). 

C. Detectability of the gravitational-wave signal 

The prospects to determine / pca k are estimated by 
computing the signal-to-noise ratio (SNR) via (SNR) 2 = 

4 J °° s^/f\ df (see e.g. [82j ]) adopting the one-sided 
spectral density Sh(/) of the strain noise of the Advanced 
LIGO detector, which has a comparable sensitivity to 
the Advanced VIRGO observatory. In Tab. [H] we list the 
SNRs for the postmerger signal and the emission during 
the whole simulation time starting from the late inspi- 
ral phase. Here we correct the underestimation of the 
GW amplitude due to the usage of the quadrupole for- 
mula (see Appendix |A| . The SNRs are computed for a 
distance of D = 20 Mpc assuming an optimal inclination 
of the binary and an optimal orientation of the detec- 
tor. Because the GW amplitude scales with 1/D, also 
the SNR is proportional to 1/D and one can easily ob- 
tain the optimal detection horizon for a chosen SNR. Re- 
quiring for instance a SNR of 2, the postmerger signal 
alone can be detected up to a distance of about 20 to 
45 Mpc depending on the EoS if we ignore results for 
EoSs which are incompatible with the observation of the 



1.97 M Q NS [3 (marked by an asterisk in Tab. HJ) and 
that yield somewhat lower SNRs. Including also the EoS 
dependent late inspiral wave train boosts the optimal de- 
tection horizon to about 150 Mpc. While a SNR of 2 is 
rather low, the preceding inspiral signal (not covered by 
our calculations) provides detailed information on im- 
portant emission features such as the time of merging, 
the distance, and the involved masses, justifying a low 
threshold SNR (in particular, the inspiral signal means a 
warning that one should search for the postmerger emis- 
sion in the close temporal vicinity). Interestingly, models 
with EoSs that result in a high SNR in the postmerger 
phase, yield in general a lower SNR for the full signal of 
the simulation and vice versa. This can be understood 
from the observation that for more compact stars the in- 
spiral phase lasts longer increasing the full SNR. In these 
cases the frequencies of the late-time GW emission are 
generally relatively high, which leads to a low SNR in 
the postmerger phase because of the reduced sensitivity 
of the GW detector at high frequencies. 

For an estimate of the detection rate with Advanced 
LIGO one considers the Milky Way Equivalent Galaxies 
(MWEGs) accessible by the detector. The reduction of 
the rate due to random source and detector orientations 
can be accounted for by dividing the optimal detection 
horizon by l/y/5 [28|]. For the chosen SNR the number 
of MWEGs equals 145 to 1190 for the postmerger sig- 
nal alone, taking advantage of the local over-density of 
potential host galaxies j28[. The inclusion of the late 
inspiral signal increases the number of probed MWEGs 
to more than 14,000. Adopting the "realistic" and the 
"high" merger rates listed in [2g|, the prospects for the 
detection of the postmerger signal alone are estimated to 
be 0.015 to 1.2 events per year. For the combined signal 
1.4 to 14 detections per year can be expected. Finally, 
we emphasize that with the increased sensitivity of the 
proposed Einstein Telescope [831 ] a significant number of 
detections with high SNRs will be possible. 

The uncertainties in the determination of the peak fre- 
quency by a GW detection can be estimated by evaluat- 
ing the Fisher information matrix for a single-parameter 
family of waveforms parametrized by / p0 ak- This ap- 
proach considers to which extent two different waveforms 
are distinguishable by a detector with a certain noise level 
(see e.g. [82|). Following the procedure of [3l[, we arrange 
our models in Tab. [TTI with increasing / pea k and compute 
for two subsequent models A and B with corresponding 
GW amplitudes h A and h B 

ffA _ fB Y2 
c r2 ^ peak J peak/ \ 

; = (h A - h B \h A - h B ) ( ) 

with the inner product between two signal hi and 

{h 1 \h 2 )=41tej — Sh(/) — ^' 

where the Fourier transforms of the signal enters (an as- 
terisk marks the complex conjugated). In Tab. [XT] the 
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entry of Sf for a given EoS corresponds to the value 
obtained by ([1]) for the waveforms of this EoS and the 
subsequent EoS in the table (e.g. for the eosL and the 
MSI EoS we compute Sf pm = 0.062 kHz). Here, a po- 
lar distance of 20 Mpc and optimal detector alignment 
are assumed. We observe uncertainties in / pca k of typi- 
cally better than 50 Hz (on average 40 Hz). Note that 
these values are smaller than the ones reported in (40| 
because a larger number of models is employed leading 
to a finer coverage of the / poa k-parameter space. It is im- 
portant to mention that Sf computed for the postmerger 
phase serves only as a rough estimate for the uncertain- 
ties in the / pea k determination because theprocedure is 
not fully applicable in the low-SNR limit |3lL HiJ . How- 
ever, we conduct the same analysis including the late 
inspiral phase of the signal, i.e. with a sufficiently high 
SNR. For a fair comparison the signal is constrained to 
/ > 1 kHz. In this case smaller uncertainties are found 
(on average 27 Hz), see column 7 in Tab. ITTT 



V. DEPENDENCE ON STELLAR 
PARAMETERS 

It is the goal of this work to establish the relations 
between features of the GW signal and properties of the 
EoS, which can be either thermodynamical properties or 
stellar parameters. To this end we determine the dom- 
inant frequency / pC ak of the postmerger GW emission 
as discussed in the previous section for all EoSs which 
lead to the formation of a DRO (sec Tab. |nj). In the 
following we investigate the dependence of / pca k on vari- 
ous NS and EoS properties, which in turn allows to con- 
strain or measure these quantities when / pea k is deter- 
mined from a GW observation. We stress that the to- 
tal binary mass and the initial mass ratio of the binary 
components can be determined from the GW inspiral sig- 
nal [201 . Throughout this and the following section we 
discuss exclusively binaries with two 1.35 M Q NSs. 



A. Neutron star radii 

In [40l ] we have identified an anticorrelation between 
the peak frequency and the minimum radius of nonspin- 
ning NSs. Fig. [5] shows this relation for our extended 
set of EoSs. All data points correspond to simulations of 
1.35 Mq-1.35 Mq binaries. The figure adopts the same 
color scheme as Fig. [4] i.e. the colors represent the four 
different possibilities for the technical implementations 
of the EoSs including the treatment of thermal effects. 
The crosses denote the results of simulations with EoSs 
consistent within the error bars with the observation of 
the 1.97 M Q pulsar [l0| . while plus signs belong to cal- 
culations employing EoSs which are excluded by this de- 
tection. The triangles correspond to the two EoSs de- 
scribing absolutely stable strange quark matter (MIT60 
in red, MIT40 in black). Note that we keep this notation 



in all figures of this paper except for Fig. [25] In Fig. [9] 
one observes a clear anticorrelation between / pea k and 
i? max with a slightly steeper slope at smaller radii. Note 
that the data points for EoSs implemented as piecewise 
polytropes (class iii) or with a high-density regime de- 
scribed by piecewise polytropes (class iv) tend to show 
some scatter from the relation established by the purely 
microphysical EoSs (black and red symbols). Since these 
models involve additional approximations or simplifica- 
tions (see Sect. Mil and reasoning below), the correspond- 
ing data points are shown with smaller symbols to em- 
phasize the correlation of the different cases. Moreover, 
excluded models with Mp ax < 1.93 M Q (lower bound of 
the error bar given by [10]) are plotted with smaller plus 
signs, although they are not extraordinary outliers in this 
figure. 

The somewhat displaced black cross (microphysical, 
barotropic EoS) at i? max = 11.48 km and / poa k=2.33 kHz 
displays the data of the Glendnh3 EoS, which has a 
strikingly different mass-radius relation in comparison 
to other microphysical EoSs (see black curve with R = 
17 km and M = 0.5 M Q in Fig. [4] or black dashed line 
in Fig. 1 of [lo]). Note that this EoS leads to a mass- 
radius relation which seems incompatible with theoretical 
calculations of the EoS of subnuclear matter [45[. Fur- 
thermore, we stress that among the purely microphysical 
models only the strange quark matter EoSs behave as 
outliers. However, if strange quark matter was the true 
ground state of matter and one of these EoSs was the cor- 
rect EoS of high-density matter, this would lead to clear 
observational signatures (e.g. in the form of nuggets of 
strange quark matter in the cosmic ray flux) which would 
discriminate this scenario from ordinary NSs composed 
of nuclconic matter (see [3, H3, El| ) . For all these reasons 
we classify the results of these three particular EoSs as 
outliers not spoiling the quality of the relation evident 
from the black and red crosses. 

Figure \W\ shows the peak frequencies from our set of 
models versus the radius Ri.s of the corresponding nonro- 
tating NS with 1.8 M . Considering in Fig. [TOlonly EoSs 
which arc in agreement with the NS mass limit of [l0( , one 
recognizes that the relation between the GW frequency 
and the NS radius is closer in comparison to Fig. [9] In 
particular, some of the EoSs modeled as piecewise poly- 
tropes (green and blue symbols) better fulfill the relation 
indicated by the purely microphysical EoSs. As will be 
discussed in Sect. IV CI the density regime probed with 
the merger remnant, i.e. affecting / poa k, corresponds 
roughly to the densities encountered in nonrotating NSs 
with masses below 1.8 Mq. 

It is striking in Fig. [10] that EoSs incompatible with 
the (1.97 ± 0.04) Mq limit (plus signs), appear as out- 
liers at smaller radii. This can be understood from the 
fact that for these EoSs the radius of a 1.8 Mq NS is 
near radius of the maximum-mass configuration because 
M max is close to 1.8 M . This means that R\.% is outside 
of the region of the mass-radius relation with approxi- 
mately constant radius and its location near the limiting 
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TABLE II: All EoS models for which the merging of two 1.35 Mq NSs results in the formation of a DRO. / pca k denotes 
the dominant GW frequency in the postmerger phase, and FWHM characterizes the corresponding width of the peak in the 
spectrum. SNRs are computed for a distance of 20 Mpc assuming an optimal source and detector orientation considering the 
postmerger signal only (SNR pm ) and the full signal above 1 kHz (SNR tot ). Uncertainties in the determination of the peak 
frequency are estimated by <5/ pm and 5/tot for the postmerger phase and the signal above 1 kHz, respectively (see main text for 
details). The last column provides the radius of a sphere enclosing a rest mass of 2.6 Mq at 8 ms after merging (see Sect. IV CI . 
EoSs marked with an asterisk are incompatible with the observation of the 1.97 Mq pulsar. 
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mass corresponds to a smaller radius (see Fig. [4j). Also 
in Fig. [TU] the data point of the Giendnh3 EoS slightly 
deviates from the anticorrclation between / pea k and i?i.s- 
In contrast to Fig. [TUJ a plot using the radius R\.& of 
a 1.6 Mq static NS as characteristic EoS property re- 
veals all EoSs fulfilling a very tight anticorrelation (see 
Fig. HTI) . Even the black cross belonging to the Glendnh3 
simulation is now located on the relation. Again, one 



finds a slight steepening of the slope at smaller radii. The 
EoSs of strange quark matter lie slightly off the relation, 
whereas in Fig. [TU] the MIT40 data point is consistent 
with the correlation (MIT60 has a relatively low M max 
and is therefore expected to appear at a smaller radius as 
argued above). In this context we stress that MIT40 and 
MIT60 describe bare strange stars. In principle, strange 
stars could carry a nuclear crust with a density below 
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the neutron drip density of 4 • 10 11 g/cm 3 (see e.g. @]). 
This crust would have only a little mass of the order of 
10~ 5 Mq and, thus, would be dynamically unimportant, 
i.e. it is unlikely to change / pea k- But the crust would 
increase the radius of a nonrotating strange star by a 
few hundred meters. This means that by considering the 
nuclear crust the triangles would move to the right, for 
instance about 500 meters for MIT40 and about 300 me- 
ters for MIT60 in the case of stars with 1.6 Mq. Hence, 
also these EoSs, in particular MIT40, are consistent with 
the relation that exists between / pca k and Ri $ for the 
nuclear EoSs. 

Finally, Fig. [T2l shows the dominant GW frequency as 
a function of the radius R1.35 of a 1.35 M© NS. In com- 
parison to Fig. [TT] the scattering is larger, although the 
relation in particular for the accepted purely microphys- 
ical models is still very tight (excluding the data point of 
the Glcndnh3 EoS). Interestingly, in Fig. [T2] outliers are 
on the opposite side of the basic relation as compared to 
their location in the plot using R max (Fig. [5]). 

Note that the different EoS classes (colors) in Figs. [9} 
[T2] sample in each case nearly the fully parameter range 
and do not indicate any trend or bias due to the different 
implementations of the EoSs aside from the increased 
scatter for implementations using piecewise polytropes. 
The EoSs of class (iv) (blue) cover a broad range of pos- 
sible behaviors at intermediate and high densities, which 
are partially very extreme (e.g., very high pressure and 
sound speed at high densities, see Figs. [5] and [3]). There- 
fore, it is expected that the resulting variations will also 
span a broad range, which is however consistent with 
the chiral effective field theory constraints at saturation 
densities and below. The models of class (iii) involve a 
twofold simplification that can explain the larger devia- 
tions from the correlations. First, the fits of the EoSs do 
not perfectly match the underlying microphysical model 
(e.g. in the sound speed, see (HI), leading to peak fre- 
quencies which may be slightly different from those ob- 
tained by the original model. Second, due to the usage 
of the fit also the stellar parameters of nonrotating NSs 
differ slightly from those obtained with the original EoSs. 
Bear in mind that the same reasoning for EoSs of class 
(iii) and class (iv) EoSs also applies to all following rela- 
tions discussed in this paper (Fig. IT3"1 to l2"Tj) . 



B. Fits and residuals 

To quantify the discussion above and to introduce a 
measure for the scatter inherent to the presented rela- 
tions, we fit power laws of the type Rtov — (/peak) b +c 
through the data points of Figs. with a, 6, c being pa- 
rameters to be obtained by a least-square fit. Rtov de- 
notes the corresponding radius of nonrotating NSs used 
as characteristic EoS parameter in Figs.l9lTT2l Because of 
the reasons mentioned in the last paragraph of Sect. IV Al 
we restrict the set of EoSs determining the parameters to 
the accepted purely microphysical models omitting also 
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FIG. 9: Peak frequency of the postmerger GW emission ver- 
sus the radius of the maximum-mass configuration of non- 
rotating NSs for different EoSs. Red symbols correspond to 
microphysical EoSs with a consistent temperature treatment 
(class i), black symbols show data points for barotropic mi- 
crophysical EoSs (without temperature and electron fraction 
dependence) (class ii) , green (smaller) symbols belong to EoSs 
implemented as piecewise polytropes fitting barotropic micro- 
physical EoSs (class iii), and blue symbols display results for 
microphysical EoSs at low densities with high-density exten- 
sions by piecewise polytropes (class iv). Classes (ii) to (iv) 
are supplemented with an ideal-gas component for mimicking 
thermal effects. Plus signs indicate EoSs which are excluded 
by the observation of a 1.97 Mq pulsar [Toj ] . EoSs describing 
absolutely stable strange quark matter are denoted by trian- 
gles. Note that the MIT60 EoS (red triangle) is ruled out by 
the 1.97 Mq mass limit. 



the EoSs of strange quark matter and the Glendnh3, 
which is in conflict with radius constraints discussed 



451 ] . The derived parameters are given in Tab. Mil for 



the four different characteristic radii of nonrotating NSs. 
Tab. Mil also quantifies the deviations of the data points 
from the corresponding fit. We list the maximum value 
and the average value of the residuals for the accepted 
models only and for all models except for the strange 
quark matter EoSs. As already apparent from Figs. [9] 
to [T2l the radius of a star with 1.6 Mq represents the 
best choice for characterizing an EoS in the sense that 
a tight relation between the dominant GW frequency of 
the postmerger remnant and the radius of the nonrotat- 
ing NS is desired. 

Considering only the 8 purely microphysical models 
(classes (i) and (ii)) that yield a peak frequency below 
2.8 kHz (89|, one finds that in Fig.[TT]the maximum resid- 
ual from a straight line is only 50 meters (average of resid- 
uals is 26 meters), allowing for an even better constraint 
of the NS radius. For the five accepted fully microphysi- 
cal models (red and black crosses) with / pca k > 2.8 kHz 
1901 we observe a maximum deviation from a linear fit of 
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FIG. 10: Peak frequency of the postmerger GW emission ver- 
sus the radius of a nonrotating NS with 1.8 Mq for different 
EoSs. Symbols have the same meaning as in Fig. [9] 



FIG. 12: Peak frequency of the postmerger GW emission ver- 
sus the radius of a nonrotating NS with 1.35 Mq for different 
EoSs. Symbols have the same meaning as in Fig. [9] 
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FIG. 11: Peak frequency of the postmerger GW emission ver- 
sus the radius of a nonrotating NS with 1.6 Mq for different 
EoSs. Symbols have the same meaning as in Fig. [9] The solid 
lines indicate the fits described by Eq. (j3}- 

only 58 meters (average of residuals is 35 meters). Hence, 
in Fig. [TT] a broken straight line represents the best ap- 
proximation to the presented scaling with 

f -0.2823 • flx.e + 6.284 for / peak < 2.8 kHz 
P ° ak ~~ \ -0.4667 • Ri 6 + 8.713 for f peak > 2.8 kHz 

(3) 

where radii arc in km and frequencies are in kHz. The 



fits are shown with solid lines in Fig.HU 

In order to quantify to which accuracy the radius of 
a 1.6 Mq NS can be obtained from a GW detection of 
the postmerger phase, we consider the maximum devia- 
tion from the derived fits and add the uncertainty in the 
measurement of the peak frequency: For / poa k we adopt 
an uncertainty of 40 Hz (see Sect. lIV Cj) . which via Eq. Q 
corresponds to 86 meters for / pea k > 2.8 kHz and 142 me- 
ters for /peak < 2.8 kHz. The larger error bar for low peak 
frequencies is caused by the flatter slope of Eq. ([3]) in this 
frequency range. We then add the maximum deviation 
of 58 meters (50 meters for / pca k < 2.8 kHz) from the 
broken line (Eq. ([3])) found for the accepted fully micro- 
physical candidate EoSs (ignoring the outliers Glcndnh3 
and MIT40 for the reasons described in Sect. IV A")) . In to- 
tal this amounts to 144 meters for / pea k > 2.8 kHz and to 
192 meters for a lower peak frequency. Being conservative 
one may add an uncertainty of up to 200 meters due to 
the use of the conformal flatness approximation subsum- 
ing also general effects by differences in the details of the 
implementation of different codes (sec Appendix |A")) . We 
do not take into account an error resulting from the ap- 
proximate treatment of thermal effects (see Appendix [B| 
because this does not apply to all EoSs. (As estimated in 
Appendix [Bl for an ideal-gas index T t h = 1.75 instead of 
2, the maximum residual of the data points from a power 
law is 126 meters.) Note that the latter two effects are 
no fundamental sources of errors and can be overcome 
by a more elaborate modeling of NS mergers and the 
high-density EoS. We also mention the possibility that 
the residuals from a fit might be further reduced by us- 
ing different radii of nonrotating NSs to characterize the 
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EoSs in different ranges of the peak frequency. 



C. Interpretation 

In a previous paper we presented arguments why the 
dominant oscillation frequency of the differentially rotat- 
ing merger remnant scales with the radius of a nonro- 
tating NS whose mass is generally smaller than the one 
of the remnant [!(| • As has been shown in [8(| , the GW 
emission at the peak frequency is generated by the funda- 
mental quadrupolar oscillation mode. The frequency of 
this mode is known to be proportional to the square root 

of the mean density, J~^js with M and R being the mass 

and the radius of the oscillating object (see d^Hfj). The 
mass of the merger remnant is approximately given by 
the total binary mass and therefore it is the same for all 
models discussed in this section neglecting small amounts 
of ejecta and differences in the inflated torus surrounding 
the central object. Hence, the peak frequency is entirely 
determined by the radius of the DRO. The radius of the 
merger remnant cannot be defined unambiguously be- 
cause one cannot identify a well defined surface of the 
object (see e.g. Fig. 5 in [37} )• Using arbitrarily the ra- 
dius of a sphere enclosing 2.6 M Q of rest mass as the 
radius of the DRO, Fig. [T3] confirms the close relation 
between / pea k and the so chosen radius of the merger 
remnant -Rremnant- Here, -Rremnant is measured 8 ms after 
merging when the oscillations of the DRO are sufficiently 
damped (see Figs. [5] and [S]). The radii of the merger rem- 
nants for different EoSs are also provided in Tab, fill (The 
data point in Fig. [13] with / pca k = 3.2 kHz, which is lo- 
cated slightly below the relation (-Rremnant = 8.53 km), 
corresponds to the fully microphysical BurgioNN EoS, 
where our somewhat arbitrary definition of -Rremnant fails. 
In particular the time, when -Rremnant is determined, is 
arbitrarily chosen. It should be sufficiently early to char- 
acterize the GW emission, but not too early when the 
DRO, and thus -Rremnant, are still strongly oscillating, 
which is the case for the BurgioNN model.) 

To understand the correlations found in Figs.Pto[T2l 
where the frequency showed a tight anticorrelation with 
radii of static TOV configurations, we hypothesize that 
for a given EoS the radius of the differentially rotating 
merger remnant of about 2.6 M scales with the ra- 
dius -Rtov of a nonrotating NS for a chosen mass. This 
hypothesis is confirmed when considering -Rremnant as a 
function of -Rtov = -R1.35, -R1.6, -R1.8, or -R max , of which 
the relation with -R1.6 shows the smallest scatter. Ne- 
glecting effects due to thermal contributions and differ- 
ential rotation, a linear relation between -R ma x and the 
radius of the most massive, uniformly rotating NS con- 
figuration was reported in [7^]. Adopting therefore a 
linear dependence between -Rremnant and .Rtov, one ex- 
pects that / pca k is proportional to R T Jy ■ When fitting a 

power law / poa k = a' • i? T ov ^° ^ ne data points of Figs. [5] 
to [T2] similar residuals as listed in Tab. IIIII are found, 
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FIG. 13: Peak frequency of the postmerger GW emission ver- 
sus the radius of a sphere enclosing 2.6 Mq of rest mass of 
the merger remnant for all fully microphysical EoSs 8 ms after 
merging. Symbols have the same meaning as in Fig. [9] 
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FIG. 14: Peak frequency of the postmerger GW emission ver- 
sus . / in geometrical units for different EoSs. Symbols 

V 1.6 

have the same meaning as in Fig. [9] 



which implies that in fact there exists a tight relation be- 
tween -Rremnant and -Rtov- Additionally, Fig. 1141 shows 
the peak frequency as a function of ./If, a quantity 

V "i.e 

which according to the above reasoning is proportional 
to the mean density of the merger remnant. This behav- 
ior is confirmed by the linear scaling evident from Fig. ll4[ 
which should be considered as an empirical finding of this 
work. 

The fact that the relation between / poa k and the radius 
of a NS with 1.6 Mq shows the best quality, can be under- 
stood by investigating the involved density regimes. For 
a given EoS the central density in the merger remnant 
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TABLE III: Parameters of power-laws fitted to the data points shown in Figs. 191121 See the main text for definitions, units are 
such that NS radii are measured in km and / pC ak in kHz. Rtov denotes which radius of nonrotating NSs is used to characterize 
the different EoSs. The last four columns provide the maximum value of the residuals and the average of the residuals for 
all accepted models of fully microphysical EoSs (without MIT60, MIT40, Glendnh3) and for all models except for the strange 
quark matter EoSs. 
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is well below the central density /3 max of the nonrotating 
maximum-mass configuration as specified in Tab. U but 
typically above the central density of the initial binary 
components, i.e., of 1.35 Mq NSs in the considered case 
(exceptions to this relatively general behavior are rather 
stiff EoSs like Shen or GS1). This can be seen by the 
density evolution in Fig. 1151 Hence, the density regime 
probed by the merger remnant and thus determining the 
GW emission is comparable to the densities encountered 
in nonrotating NSs with masses above 1.35 Mq but be- 
low M max . On the one hand this means that the value 
of i? max as a characteristic EoS property is determined 
by a density regime of the EoS that does not occur in 
the merger remnant. This explains the larger scatter in 
Fig. |U On the other hand, the radius -R1.35 depends 
only on the EoS at relatively low densities, whereas the 
structure of a DRO of a merger with M tot — 2.7 Mq is 
determined by the EoS at higher densities. Thus, the 
relatively large deviations in Fig. [12] are a consequence of 
the inability of a 1.35 Mq NS to capture the EoS behav- 
ior at densities realized in the merger remnant. 



D. Additional dependences 

Though we identify NS radii as the crucial character- 
istic EoS-dependent quantity that affects the late-time 
GW emission, additional constraints on NS properties 
by GW signal features can be derived from our survey 
of NS mergers for a large set of EoSs. Fig. \W\ shows 
the peak frequency as a function of the maximum mass 
Afmax of nonrotating NSs. The vertical line indicates the 
mass limit set by the observation of the 1.97 Mq pul- 
sar [lfj. Note that some EoS on the left of the vertical 
line (crosses) are not considered as excluded because they 
are compatible with the mass measurement within the er- 
ror bars. Figure [TBI shows that high peak frequencies tend 
to exclude high values of M max . Note that the three blue 
symbols with very high M max belong to EoSs of class (iv) 
(with piecewise-polytrope extension at high densities), 
and that these three particular EoSs are subject to an ex- 
treme stiffening of the EoS at high densities (see Figs. [2] 
and [3]). Considering our set of representative EoSs one 
can conclude that for a given peak frequency / poa k the 
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FIG. 15: Evolution of the maximum rest-mass density for the 
models with the Shen EoS (black), the Sly4 (red) EoS and 
the MIT60 EoS (blue). At the beginning of the evolution the 
shown maximum rest-mass density corresponds to the central 
density in the inspiraling NSs, while after merging the maxi- 
mum value of the rest-mass density is found in the center of 
the DRO. 



maximum mass of nonrotating NSs can be constrained 
to roughly M max < 4.25 Mq - O.5(M /kHz) ■ / pcak (see 
dashed line in [TO)) . 

In addition to the peak frequency one may also con- 
sider to use the width of the peak to characterize the 
GW emission during the postmerger phase, although it 
is unlikely to be obtained precisely from a GW observa- 
tion by Advanced LIGO or Advanced Virgo. However, 
the situation may change with the increased sensitivity 
of the Einstein Telescope. Since the peak frequency is de- 
termined by the size of the merger remnant, the width of 
the peak can be considered as a measure for the change of 
the structure of the DRO and its evolution. In the follow- 
ing we quantify the width of the peak by the FWHM (see 
Tab. [EI]). Fig. [T7I illustrates the dependence of the FWHM 
on Mmax, whereas in Fig.[TH]thc peak width is shown as a 
function of the radius of a nonrotating NS with 1.6 Mq. 
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FIG. 16: Peak frequency of the postmerger GW emission ver- 
sus the maximum mass M max of nonrotating NSs for different 
EoSs. Symbols have the same meaning as in Fig. [9] The ver- 
tical solid line marks the 1.97 M® mass limit set by [lj}. The 
dashed line represents an upper limit of M max for fully micro- 
physical EoSs (see main text). 



Using i? max to characterize the different models results 
in a picture qualitatively similar to Fig. 1181 It is striking 
that the FWHM does not show a clear dependence on the 
stellar parameters characterizing a given EoS, instead for 
a given FWHM a large spread in M max and Ri,e is found. 
A slight trend is noticeable that EoSs with high M max or 
larger Ri.q tend to have smaller peak widths. Further- 
more, one recognizes that only models that are excluded 
by the observation of the 1.97 Mq pulsar yield increased 
FWHMs, though even in these cases the absolute values 
are not higher than about 250 Hz (except for the strange 
quark matter MIT60 EoS). However, we emphasize that 
not all excluded models show this behavior. 

The broadest peak with FWHM=460 Hz is found 
for the MIT60 EoS describing absolutely stable strange 
quark matter. It seems unlikely that this finding is a 
particular consequence of the fact that this EoS models 
strange quark stars because for MIT40 a narrow peak is 
found. Rather, what distinguishes MIT60 from the other 
EoSs is the early occurrence of the delayed BH forma- 
tion taking place about 3 ms after merging, i.e. during a 
period when the postmerger remnant is still vividly emit- 
ting GWs. This leads to the conclusion that the broad- 
ening of the peak is connected to the dramatic structural 
changes in the DRO prior to collapse, visible for example 
in the evolution of the central density (see Fig. [TBI . Sec 
also the discussion in Sect. IVIII 

We stress that MIT60 is the only model out of the 
38 discussed here, where the delayed collapse occurs in 
the period of strong GW emission, i.e. at a time when 
the peak in the GW spectrum is shaped. In the re- 



maining simulations either the prompt collapse to a BH 
takes place or the forming DROs are stable for at least 
10 to 15 ms until the GW emission has significantly de- 
creased. Only in one other case, namely the H3 EoS, we 
observe the delayed collapse before the end of the sim- 
ulation about 25 ms after merging. Note that MIT60 
yields an M max by which this EoS is excluded because 
of M max < 1.97 Mq. The broadest peak among the ac- 
cepted models of our survey is only 205 Hz (for the LS220 
EoS). 

Considering these observations and the large variety in 
our representative sample of EoSs we conclude that the 
early occurrence of a delayed collapse of the DRO is very 
unlikely to happen and, therefore, a broadening of the 
peak in the postmerger GW spectrum seems improbable 
to occur at least for the binary setup considered here. 
Beside this empirical finding, to our knowledge there is 
no case reported in the literature where the merger of two 
stars with 1.35 Mq yields an extraordinary broad peak 
in the GW spectrum for an EoS which is compatible with 
the NS mass measurement of [Toj . 

In [38| the delayed collapse within the first 5 ms of the 
lifetime of the DRO was found to be generic in a certain 
range of total binary masses slightly below the threshold 
for the direct BH formation. In view of this it is unex- 
pected that we do not observe any case of such an early 
delayed collapse except for the excluded strange quark 
matter MIT60 EoS when we consider our very large sam- 
ple of EoSs. For instance, according to [38[ the merger of 
two NSs of 1.35 M with the Sly4 EoS should have led 
to a collapse within 5 ms after merging. To explore this 
difference we perform additional calculations for the Sly4 
EoS with slightly increased binary masses. The result- 
ing lifetimes of the remnants are summarized in Tab. IIVI 
for simulations using T t h — 2 (our standard choice; see 
Sect. |TT| and T t h = 1-357 as employed in [38|. 

For r t h = 2 we find either the prompt collapse to a BH 
or the formation of relatively long-lived merger remnants 
depending on the total binary mass. In contrast, with a 
description of thermal effects that corresponds to a very 
inefficient shock heating (r t h = 1.357) we indeed identify 
a parameter range for which short-lived merger remnants 
with lifetimes below 5 ms occur. This is in qualitative 
agreement with the findings of [38}. Our calculations 
thus suggest that the frequent occurrence of short DRO 
lifetimes (below 5 ms) in [38| is a consequence of the 
assumption of a very low shock-heating efficiency and it 
shows that the choice of r t h has a crucial influence on the 
stability of the merger remnant. This is understandable 
since thermal pressure support has generally a stabilizing 
effect. (In [3 the thermal ideal-gas index was extracted 
from microphysical EoSs including a consistent descrip- 
tion of temperature effects, showing that T t h in the range 
from 1.5 to 2 represents a suitable approximation. Sec 
Fig. 2 in [13.) 

The quantitative differences in Tab. IIVI might be ex- 
plained by differences in the exact implementation of the 
EoS (original table vs. piecewise polytropic approxima- 
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FIG. 17: FWHM of the peak structure of the postmerger GW 
emission versus the maximum mass M max of nonrotating NSs 
for different EoSs. Symbols have the same meaning as in 
Fig.El 
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FIG. 18: FWHM of the peak structure of the postmerger GW 
emission versus the radius of nonrotating NSs with 1.6 M@ for 
different EoSs. Symbols have the same meaning as in Fig. [9] 



tion of the EoS) and the different degree of sophistication 
of the simulations (SPH in combination with the con- 
formal flatness approximation and a GW backreaction 
scheme vs. grid-based hydrodynamics in full general rel- 
ativity where the simulations start more than five orbits 
before merging) . See also Appendix [A] 

In the context of DRO lifetimes we also mention that a 
slight decrease in the total binary mass (M tot — 2.68 M© 
instead of M to t = 2.7 M Q ) for the model with the eosAU 
prevents the prompt collapse and instead leads to the for- 
mation of a DRO which is stable for at least 14 ms and 
generates a narrow peak of 169 Hz at / pca k = 3.94 kHz. 
Hence, based on this marginal case and on our survey 
we conclude that the generic outcome of a NS merger 
of two stars with 1.35 Mq is either the prompt collapse 
or a relatively long-lived DRO (with a lifetime of typi- 
cally more than 5 ms) where shock heating is crucial for 
the stability of the DRO. We stress that for a thorough 
investigation of the stability and lifetime of merger rem- 
nants the consideration of a large sample of temperature 
dependent EoSs is essential. 



VI. DEPENDENCE ON THERMODYNAMIC AL 
PROPERTIES 

Section [V] provides evidence that NS radii play a cru- 
cial role for determining the GW emission of NS mergers. 
Stellar parameters can be considered to be a consequence 
of the bulk properties of the high-density EoS because 
they are integral values associated with the solution of 
the TOV equations taking into account the EoS behavior 
over a large density regime. In this section we directly re- 
late the dominant GW frequency of the postmerger phase 



to specific properties of the NS EoS (see Tab. HJ. Fig. [19] 
shows a clear correlation between / pea k and the central 
energy density of the nonrotating maximum-mass config- 
uration of the TOV solutions. However, in comparison 
to Figs. |H]to[T2]onc finds a wider scatter and one cannot 
infer a clear functional dependence. At most exclusion 
regions can be identified. As argued in Sect. [V] the devi- 
ations are understandable because only a certain density 
range, which does not include e max , is probed by the 
postmerger remnant. This effect is most pronounced for 
hyperon EoSs, where in contrast to nucleonic matter the 
behavior of the high-density EoS is changed by the ap- 
pearance of hyperons. Consequently, the hyperon EoSs 
are found on the lower right of the band established by 
the correlation of other EoSs. 

Figure [201 which displays / pea k as a function of the 
pressure at 1.85po confirms these arguments. Here, the 
intrinsic scatter is smaller because one considers a charac- 
teristic quantity at a moderate density which is realized 
in the merger remnants. One notices that a low pres- 
sure at the fiducial rest-mass density of 1.85po causes a 
higher peak frequency. This is plausible because a re- 
duced pressure yields more compact stellar objects re- 
sulting in higher oscillation frequencies as discussed in 
Sect.[V] A direct correlation between P(1.85/?o) and the 
radius of nonrotating NSs has been pointed out in [l| and 
is confirmed by the values given in Tab. [I] 

Using the speed of sound as a characteristic EoS prop- 
erty reveals a relation as shown in Fig. [5T] The sound 
speed is determined at 1.85po- At high peak frequencies 
above 3 kHz the anticorrclation is rather tight, whereas 
the chosen EoS feature cannot be constrained signifi- 
cantly in cases of low peak frequencies. Note that using 
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TABLE IV: Lifetime in ms of the postmerger remnant in simulations for the Sly4 EoS. In the first row the masses of the 
simulated binary components are given in Mq . The entry "prompt BH" indicates direct BH formation. Bars mark cases where 
no calculations have been performed. Note that we report lower bounds on the lifetime when the collapse to a BH does not 
occur until the end of the simulation. 
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FIG. 19: Peak frequency of the postmerger GW emission ver- 
sus the maximum central energy density of stable nonrotating 
NSs for different EoSs. Symbols have the same meaning as in 
Fig. El 
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FIG. 20: Peak frequency of the postmerger GW emission ver- 
sus the pressure at 1.85 times nuclear saturation density for 
different EoSs. Symbols have the same meaning as in Fig. [9] 



the speed of sound at saturation density reduces the scat- 
ter at frequencies below 2.8 kHz but leads to an increase 
of the variations at high / pea k (not shown). As a general 
finding we stress that choosing other fiducial densities for 
the pressure or the sound speed does not yield tighter cor- 
relations over the whole parameter range. We also note 
that the locus of the strange quark matter EoSs in dia- 
grams like Figs. [20] or [21] is highly sensitive to the chosen 
reference density (see Figs. [5] and [3] for the determination 
of the characteristic properties depending on the fiducial 
density). 



VII. VARIATION OF BINARY PARAMETERS 



A. Dependence on stellar radii 
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In the previous sections NS mergers with a total binary 
mass of 2.7 M to t and a mass ratio q = M1/M2 = 1 were 
considered and very tight anticorrelations in particular 
between the dominant GW frequency in the postmerger 
phase and the stellar radii of nonrotating NSs were ob- 



FIG. 21: Peak frequency of the postmerger GW emission ver- 
sus the sound speed at 1.85 times nuclear saturation density 
for different EoSs. Symbols have the same meaning as in 
Fig.H 
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served. In the present section we explore these relations 
for different binary setups. In our discussion we include 
the results of simulations of binaries with total system 
masses (i.e. the gravitational mass of the two compo- 
nents in isolation added) of 2.4 Mq, 2.7 M and 3.0 M 
each with a mass ratio 5=1. Calculations arc performed 
for all fully microphysical EoSs (classes i and ii). For the 
EoSs with a consistent temperature treatment (class i) 
we also test the effect of asymmetric initial masses by 
simulating 1.2-1.5 Mq mergers. All results are given in 
Tabs. [V] (symmetric binaries) and IVII (asymmetric bina- 
ries). Figs. [2H US] and [M] display the peak frequency of 
the postmerger stage as a function of the radii of nonro- 
tating NSs with 1.35 M , 1.6 M and 1.8 M for the 
aforementioned binary setups analogous to Figs.fT0lto[T2l 
The simulations with M tot = 2.7 Mq, which were al- 
ready discussed above, are shown with the same symbols 
used previously, i.e. crosses and plus signs. The same 
color scheme as in the previous sections is adopted. Low- 
mass binary merger are plotted with squares. The results 
of calculations with M to t = 3.0 M are given by dia- 
monds. Data points corresponding to excluded EoSs are 
displayed by smaller symbols. Simulations with q ^ 1 arc 
indicated by circles. Note that the strange quark-matter 
cases are highlighted by special symbols only for the 1.35- 
1.35 Mq mergers. However, the symbols belonging to a 
certain EoS can be readily identified by their location at 
the same NS radius. 

The M to t = 2.7 M simulations were extensively dis- 
cussed in Sect. [V] Considering the calculations with two 
1.2 Mq stars one finds in Figs . |2"21 to |2"41 that also for this 
binary setup relations between / poa k and the NS radii ex- 
ist. The anticorrelation found with the radius -R1.35 is the 
tightest in comparison to the others. This is fully consis- 
tent with the interpretation in Sect. IV Cl that the merger 
remnant probes mostly a certain density regime which 
for lower Af to t corresponds to the radii of lower-mass 
NSs. Consequently, the scatter is enhanced when using 
the radii of nonrotating high-mass NSs as characteristic 
properties (Figs. [23l and [24]) . The strange quark matter 
EoSs are off the relations in all three plots. In addition, 
only the result for the Glendnh3 EoS (#1.35 = 14.52 km) 
seems to deviate slightly from the very tight anticorrela- 
tion established by the other EoSs in Fig. [22] 

Choosing Ri.g as characterizing EoS property in 
Fig. [231 all excluded models (small symbols) occur as 
outliers at smaller radii but seem to fulfill their own re- 
lation. This behavior is in analogy to the one in Fig. [10] 
(see discussion in Sect. [V]) and can be explained by the 
fact that for these EoSs the NSs with 1.6 M are not lo- 
cated on the nearly vertical sequence of the mass-radius 
relation as it is the case for EoSs with higher M max . 

In Fig. [23] the two EoSs appearing at nearly the same 
radius of about 13 km are the SKA EoS and eosO (black 
symbols). While they show some difference in / pca k for 
the M tot = 2.4 M Q and M tot = 3.0 M runs, the 
peak frequencies of the 1.35-1.35 M mergers are very 
close. This illustrates the importance of choosing an 



appropriate characterizing NS radius dependent on the 
total binary mass, cf. the location of these EoSs in 
Figs. [22] and [24j Accordingly, for mergers of two stars 
with 1.5 Mq the radius of a nonrotating 1.8 Mq NS 
leads to the most narrow anticorrelation, although com- 
pared to the simulations with M to t =2.4 M or 2.7 M 
still larger scatter is visible. However, note that only the 
three models based on barotropic EoSs (black diamonds) 
seem to yield peak frequencies that are systematically on 
the lower side while the remaining models form a clear 
relation (that can be described for instance by a power 
law). The behavior of these three class- ii EoSs at smaller 
radii may be an artifact connected to the approximate 
temperature treatment needed for these models. The us- 
age of r t h = 2 may underestimate / pea k hi comparison 
to a consistent description of thermal effects (sec Ap- 
pendix [B] and (l9j ) , which becomes in particular impor- 
tant for high-mass binaries and EoSs yielding compact 
NSs. 



Exploring the influence of M to t on the peak frequency 
for different EoSs, a stronger effect is found for mod- 
els with compact NSs, i.e. at higher GW frequencies. 
While for cosL, for example, / pea k ranges from 1.76 kHz 
to 1.98 kHz when M tot varies from 2.4 M to 3.0 Mq, 
differences in the peak frequency exceeding 500 Hz are 
observed for EoSs with smaller NS radii for the same 
variations in the total binary mass. In comparison, the 
effect of the initial mass ratio of the binaries is more mod- 
est. The peak frequencies of 1.2-1.5 M mergers are in 
general slightly lower than the corresponding / pea k of the 
symmetric setup with the same total binary mass. This 
trend can be expected because for an asymmetric binary 
during the coalescence more matter is spread into the 
torus surrounding the central object. Hence, the DRO is 
less massive and thus oscillates at lower frequencies. The 
increased frequency for the 1.2-1.5 M merger with the 
excluded LSI 80 EoS is a consequence of the gravitational 
collapse only 9 ms after merging, which is accompanied 
by a compactification of the remnant prior to the col- 
lapse. Note that in Figs. [22]to[M]wc d° not include the 
result of the asymmetric merger with the MIT60 EoS be- 
cause this model collapses about 2 ms after merging. For 
the accepted models the deviations of / pea k from the re- 
sults of symmetric binaries seem to increase for smaller 
NS radii, but the GS1 EoS makes an exception to this 
tendency. The larger differences for smaller radii can be 
understood from the finding above that changes in the 
mass of the central DRO have a stronger effect on / pea k 
for such models. In summary, the influence of the ini- 
tial mass ratio is relatively small, at most 90 Hz for the 
accepted models of class i. This demonstrates the insen- 
sitivity of the presented relations to small variations of 
the initial mass ratio q and shows the robustness of our 
method for constraining NS radii with respect to uncer- 
tainties in the accurate determination of q from the GW 
inspiral signal. 
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FIG. 22: Peak frequency of the postmerger GW emission ver- 
sus the radius of a nonrotating NS with 1.35 Mq for differ- 
ent EoSs. The same color scheme as in Fig. [9] is adopted. 
Squares mark the results of simulations with M to t = 2.4 Mq, 
crosses and plus signs indicate the peak frequencies of merg- 
ers of two stars with 1.35 Mq, circles correspond to calcula- 
tions of 1.2-1.5 Mq mergers, and diamonds denote the out- 
come of symmetric binary coalescences with M to t = 3.0 Mq. 
Triangles belong to the 1.35-1.35 Mq simulations with the 
strange quark matter EoSs (MIT40 and MIT60). Note that 
the other binary setups for these EoSs are not specially high- 
lighted. Small symbols indicate that the corresponding EoS 
is excluded by the 1.97 Mq pulsar observation (except for the 
strange quark matter EoS MIT40, which is not excluded). 
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FIG. 23: Peak frequency of the postmerger GW emission ver- 
sus the radius of a nonrotating NS with 1.6 Mq for different 
EoSs and varied binary setups. Symbols have the same mean- 
ing as in Fig. [22] 



FIG. 24: Peak frequency of the postmerger GW emission ver- 
sus the radius of a nonrotating NS with 1.8 Mq for different 
EoSs and varied binary setups. Symbols have the same mean- 
ing as in Fig. [22] 



B. Prompt collapse 

It can be noticed in Figs.[22]to[24]that in comparison to 
the 1.35-1.35 A/0 simulations there are more data points 
for 1.2-1.2 AIq binaries and less for 1.5-1.5 Mq mergers 
because with higher total binary masses for more EoSs 
the prompt collapse to a BH occurs. Characterizing a 
given EoS by M max and -R1.35, Fig. [231 illustrates the sit- 
uation: Models are marked with asterisks there when the 
direct BH formation is observed for 1.35-1.35 Mq merg- 
ers. Crosses indicate EoSs which lead to the prompt col- 
lapse for a total binary mass of M to t =3.0 Mq. The re- 
maining fully microphysical models of 1.5-1.5 Mq merg- 
ers with DRO formation are denoted by dots. Note that 
we do not find cases where the coalescence of two NSs 
with 1.2 Mq results in the direct occurrence of a BH. The 



data point at 



10.43 km and M n 



1. 



Mq 



marks the MIT60 EoS, which as a strange quark matter 
EoS should be disregarded in this discussion. Despite the 
coarse sampling in the -ffi.35 — M max space one can iden- 
tify the region in the parameter space that allows for the 
direct BH formation for a given M tot . For a fixed total 
binary mass the boundary between this scenario and the 
formation of a DRO may be roughly approximated by 



M max = -0.6(^)^.35 
km 



M c 



l [so t 



(M tot ), (4) 



where only the offset depends on M to t- We determine 
M offset (2.7 Mq) = 8.45 M Q and M offsct (3.0 Mq) = 



9.3 Mq (see dashed lines in Fig. 



These findings 



are qualitatively consistent with the analysis of (381 ] , who 
proposed for the threshold binary mass distinguishing the 
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FIG. 25: Equation of state properties for direct BH forma- 
tion are characterized by the corresponding maximum mass 
of nonrotating NS and the radius of a nonrotating NS with 
1.35 Mq. Models for which the merger of a symmetric binary 
with A/tot =2.7 Mq leads to the prompt formation of a BH 
are displayed with asterisks. Crosses indicate cases where the 
coalescence of two stars with 1.5 Mq each results in a prompt 
collapse. Dots denote simulations where the formation of a 
DRO is found for binaries with M to t = 3.0 Mq. Note that we 
only consider EoSs of classes i and ii. The small cross marks 
the MIT60 EoS describing absolutely stable strange quark 
matter. Models collapsing promptly for Mtot = 2.7 Mq also 
form directly a BH in a 1.5-1.5 Mq merger. The dashed lines 
indicate the approximate boundary between prompt and de- 
layed collapse for M to t = 2.7 Mq and M to t = 3.0 M© (see 
text). 



two scenarios of a prompt and a delayed collapse a depen- 
dence Mthreshoid = kM max . The parameter k was found 
to depend on the radius of a 1.4 Mq NS, which is con- 
firmed by our calculations. A quantitative comparison 
is not meaningful because of the coarse sampling of our 
models in M to t and because of the usage of r t h = 1.357 
in [lU which implies a very low pressure increase associ- 
ated with shock heating. But we stress that the relation 
between fc and R\a is not very tight, which can be also 
seen in the data of 



This impedes the determination 
of M max by means of this scaling. Besides, there are in- 
dications for k being also sensitive to M to t- Note that in 
this subsection we focus entirely on symmetric mergers 
(q = 1). The picture may change for asymmetric setups. 



the peak in the GW spectrum. This is demonstrated in 
Fig. [25] showing the FWHM as a function of the radius 
iJi.6 of a NS with 1.6 Mq adopting the same notation as 
in Fig. [221 As already indicated in Fig.[TH]the simulations 
with EoSs yielding more compact NSs result in broader 
peaks, which is confirmed by Fig. [25] also for other total 
binary masses. It is apparent that in calculations with 
excluded EoSs (small symbols) generally larger FWHMs 
are observed. Asymmetric binary setups also tend to 
generate broader peaks (circles), which for many EoSs 
is a stronger effect than varying M tot by 0.3 Mq. Note 
that extraordinarily broadened peaks occur only in sim- 
ulations where the DRO collapses in the phase that still 
produces a substantial emission of GWs. This is the case 
for the models LS180-1.2-1.5, MIT60-1. 2-1.5, LS180-1.4- 
1.4, and LS220-1.5-1.5. From these general observations 
we conclude that the width of the peak in general is a 
measure of the proximity of the DRO to collapse, prior 
to which strong structural changes of the remnant occur. 
This can be seen, for instance, by the outcome and the 
FWHM for all simulations with the LS180 EoS when we 
include in addition the calculation for a symmetric bi- 
nary with Mtot = 2.8 Mq (not considered in the plots 
but listed in Tab. W\ . This also means that our study 
cannot confirm that a broadened peak is a particular and 
generic consequence of the appearance of hyperons as it 
has been concluded in [3^] based on the investigation of 
only two EoSs. These differ not only with respect to their 
hypcron content but in particular also in the resulting 
stellar parameters, both of which can have crucial influ- 
ence on the dynamical behavior. For instance we find 
broad peaks also for purely nuclconic EoSs which lead 
to a collapse a few ms after merging, while not every 
EoS including hyperons causes an enlarged width of the 
peak (cf. models Glendnh3-1.2-1.2, Glendnh3-1.35-1.35, 
H3-l.35-l.35, H4-1. 35-1.35, Glendnh3-1.5-1.5). We stress 
that for the considered binary configurations broadened 
peaks are rather rare, in particular with EoSs compatible 
with the observation of a 1.97 Mq pulsar. Hence, it is 
unlikely that the broadening of the dominant peak in the 
postmerger GW spectrum jeopardizes the determination 
of NS properties by the relations presented in this work, 
because wide peaks occur only in the seldom cases when 
the gravitational collapse of the DRO sets in during the 
phase of strong GW emission. Moreover, for the models 
examined here the peak frequency of broadened peaks 
still fulfills the discussed relations. 



VIII. CONCLUSIONS 



C. Width of the peak 

Tables [V] and I VII list the results of all simulations dis- 
cussed in this section. In particular, the tables provide 
also the widths of the peaks. As a general trend one finds 
that the higher the total binary mass is, the broader is 



In this work we have investigated an extended set of 
NS merger simulations considering in total 38 microphys- 
ical EoSs. This representative sample of high-density 
EoS candidates shows a large diversity in its properties, 
e.g. concerning the resulting stellar parameters, where 
the maximum NS mass varies from 1.79 to 3.00 Mq 
with corresponding radii ranging from 8.65 to 14.30 km 
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TABLE V: Models with varied total binary masses and an ini- 
tial mass ratio of q = 1. Mi and M2 refer to the masses of the 
binary components. In the second column the ideal-gas index 
for mimicking thermal effects is given, "full" denotes a consis- 
tent temperature treatment. The third column lists the dom- 
inant GW frequency / pea k in the postmerger phase, "prompt 
BH" indicates the prompt BH formation for the given model. 
The width of the peak is provided in the last column. EoSs 
which are excluded by the pulsar mass measurement of [T3| 
are marked by an asterisk. 



EoS Mi - M 2 


Tth 


/peak [kHz] 


FWHM [kHz] 


Sly4 1.2-1.2 


2 


2.94 


0.139 


APR 1.2-1.2 


2 


3.15 


0.143 


FPS 1.2-1.2 * 


2 


3.41 


0.186 


BBB2 1.2-1.2 * 


2 


3.15 


0.183 


Glendnh3 1.2-1.2 


2 


2.12 


0.185 


eosAU 1.2-1.2 


2 


3.51 


0.244 


eosC 1.2-1.2 * 


2 


2.81 


0.196 


eosL 1.2-1.2 


2 


1.76 


0.124 


eosO 1.2-1.2 


2 


2.48 


0.111 


eosUU 1.2-1.2 


2 


3.21 


0.127 


eosWS 1.2-1.2 * 


2 


3.35 


0.170 


SKA 1.2-1.2 


2 


2.40 


0.129 


Shen 1.2-1.2 


full 


2.02 


0.143 


LS180 1.2-1.2 * 


full 


2.77 


0.195 


LS220 1.2-1.2 


full 


2.55 


0.166 


LS375 1.2-1.2 


full 


2.24 


0.130 


GS1 1.2-1.2 


full 


1.97 


0.129 


GS2 1.2-1.2 


full 


2.34 


0.126 


MIT60 1.2-1.2 * 


full 


3.01 


0.207 


BurgioNN 1.2-1.2 


2 


2.85 


0.157 


MIT40 1.2-1.2 


1.34 


2.35 


0.161 


LS180 1.4-1.4 * 


full 


3.87 a 


0.475 


Sly4 1.5-1.5 


2 


prompt BH 




APR 1.5-1.5 


2 


prompt BH 




FPS 1.5-1.5 * 


2 


prompt BH 




BBB2 1.5-1.5 * 


2 


prompt BH 




Glendnh3 1.5-1.5 


2 


2.64 


0.217 


eosAU 1.5-1.5 


2 


prompt BH 




eosC 1.5-1.5 * 


2 


prompt BH 




eosL 1.5-1.5 


2 


1.98 


0.137 


pnsO 1^1^ 

CUSW J..U - ±.0 


9 


z.ou 


n 1 rr 


eosUU 1.5-1.5 


2 


prompt BH 




eosWS 1.5-1.5 * 


2 


prompt BH 




SKA 1.5-1.5 


2 


3.01 


0.200 


Shen 1.5-1.5 


full 


2.39 


0.160 


LS180 1.5-1.5 * 


full 


prompt BH 




LS220 1.5-1.5 


full 


3.50 6 


0.326 


LS375 1.5-1.5 


full 


2.57 


0.142 


GS1 1.5-1.5 


full 


2.27 


0.130 


GS2 1.5-1.5 


full 


2.83 


0.176 


MIT60 1.5-1.5 * 


full 


prompt BH 




BurgioNN 1.5-1.5 


2 


prompt BH 




MIT40 1.5-1.5 


1.34 


2.83 


0.280 



0.6r 
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FIG. 26: FWHM of the peak structure of the postmerger GW 
emission versus the radius of nonrotating NSs with 1.6 Mq 
for different EoSs and varied binary setups. Symbols have the 
same meaning as in Fig. [22] 



TABLE VI: Models with varied initial mass ratio and M to t = 
2.7 Mq. Mi and M% refer to the masses of the binary compo- 
nents. In the second column the ideal-gas index for mimicking 
thermal effects is given, "full" denotes a consistent temper- 
ature treatment. The third column lists the dominant GW 
frequency / pea k in the postmerger phase. The width of the 
peak is provided in the last column. EoSs which are excluded 
by the pulsar mass measurement of [l^] are marked by an 
asterisk. 



Shen 1.2-1.5 


full 


2.19 


0.169 


LS180 1.2-1.5 * 


full 


~3.4 a 


~0.3 


LS220 1.2-1.5 


full 


2.80 


0.260 


LS375 1.2-1.5 


full 


2.38 


0.120 


GS1 1.2-1.5 


full 


2.03 


0.190 


GS2 1.2-1.5 


full 


2.49 


0.239 


MIT60 1.2-1.5 * 


full 


~3A b 


~0.6 



"collapse after 9 ms 
''collapse after 2 ms 



(Fig. 01). We emphasize that no selection procedure is 
applied to the choice of the EoSs (except for requiring 
M max > 1-97 Mq because of the pulsar mass measure- 
ment of Our models also include two EoSs describ- 
ing absolutely stable strange quark matter. 

For 34 EoSs of our sample the coalescence of two stars 
with 1.35 M & leads to the formation of a DRO, whose os- 
cillations generate GW emission mainly in the kHz range. 
In particular, the fundamental quadrupolar fluid mode 
produces a pronounced peak between 1.84 and 3.73 kHz 
depending on the employed EoS. Except for one strange 



"collapse after 7 ms 
'collapse after 9 ms 
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quark matter EoS, which is incompatible with observa- 
tions, the peaks are very narrow with FWHMs of typ- 
ically below 200 Hz. This pronounced feature of the 
postmerger GW emission is measurable up to a distance 
of 20 to 45 Mpc with a SNR of 2. Including also the 
late inspiral signal shifts the detection horizon up to 
150 Mpc. Assuming a reasonably optimistic merger rate 
(1000 MWEG" 1 Myr" 1 ) this implies a detection rate for 
Advanced LIGO of about 1 event per year (up to 14 with 
the inclusion of the late inspiral phase) . 

The frequency of the peak in the GW spectrum turns 
out to be a direct measure for the size of the merger 
remnant. Furthermore, the peak frequency is found to 
anticorrelate very well with the radii of nonrotating NSs 
(Figs. l9"ltoll2p. The relation is in particular tight when 
choosing the radius of a 1.6 M Q NS as a property char- 
acterizing an EoS (Fig. ITT]) ■ For such a configuration the 
density regimes encountered in the massive DRO and in 
the less massive, nonrotating NS are similar. In this case 
the maximum offset from a broken linear fit describing 
the relation between the peak frequency and the radius is 
60 meters for all fully microphysical EoSs of our sample 
which arc in agreement with current knowledge of NSs 
(mass limit set by the observation of the 1.97 M Q pul- 
sar). Considering in addition that the peak frequency 
can be measured with an uncertainty of about 40 Hz, NS 
radii can be determined by a GW detection with an accu- 
racy of about 144-200 meters depending on the peak fre- 
quency. This should be compared with the 1 km precision 
found for the radius determination from the GW inspiral 
signal with an optimal detection horizon of 100 Mpc j3lj . 
Our method as presented in this paper therefore promises 
more accurate results but with a lower detection rate. In 
addition, the robustness concerning small uncertainties 
in the initial mass ratio of the binary has been shown 
(Figs. [22] to [24]). It is also worth mentioning that inspiral 
signals will provide NS radius information of the inspi- 
raling stars, thus typically for NSs with 1.35 Mq, which 
are expected to be the most abundant cases in NS bi- 
naries. Using the postmerger signal of the same binary 
will constrain the radius of a more massive star (-Ri.6 f° r 
Mtot = 2.7 Mq), which seems advantageous because it 
probes higher densities of the supranuclear EoS. 

Examining the correlations between the dominant GW 
frequency of the postmerger signal and NS radii, we find 
that EoSs relying on piecewise polytropcs show a larger 
scatter in these relations, which is probably a conse- 
quence of the involved simplifications. In any case it is 
assuring that also with these models the scatter in the 
observed relations is still small (maximum deviation of 
357 meters). 

Concerning the strange quark matter models we note 
that they approximately fulfill the relations established 
for the nucleonic EoSs, in particular, when including 
an additional nuclear crust of strange stars. However, 
the scenario of absolutely stable strange quark matter 
is observationally distinguishable from ordinary neutron 
stars, e.g. by strangclcts in the cosmic ray flux. Whether 



strange stars actually obey their "own" peak frequency- 
radius relation should be explored if more evidence arose 
for the existence of these objects. 

As an advantage compared to the analysis of the GW 
inspiral signal, the GW peak frequencies extracted from 
our simulations can also be set into relation to other NS 
and EoS properties. By this, exclusion regions for the 
maximum mass of nonrotating NSs are identified for cer- 
tain values of the peak frequency. High frequencies tend 
to exclude high maximum NS masses and a single de- 
termination of the peak frequency may significantly con- 
strain the maximum NS mass (Fig. [16]) . This provides an 
alternative to estimating the maximum mass via the de- 
termination of the binary mass at the boundary between 
prompt and delayed collapse, which requires several de- 
tections and an a priori knowledge of NS radii [13, HH . 

Moreover, we found that the dominant GW frequency 
of the DRO scales with the maximum central energy den- 
sity of nonrotating NSs and with the pressure and the 
speed of sound at 1.85 times the nuclear saturation den- 
sity (Figs. [T9"l to l2"Tj) . All these dependences can be used 
to constrain the high-density EoS, although in compari- 
son to stellar radii as consequence of the bulk properties 
of the high-density EoS the relations with local EoS prop- 
erties are less tight. Finally, we note that the widths of 
the peaks in the GW spectra show the general trend of 
being smaller in simulations with EoSs which yield higher 
maximum masses or less compact NSs (Figs. [T71 and ITS"]) . 

Variations with the binary parameters were investi- 
gated and scalings between the peak frequency and NS 
radii were determined for simulations with total binary 
masses of 2.4, 2.7 and 3.0 Mq. Bear in mind that the 
binary masses are measurable from the GW inspiral sig- 
nal. It was noticed that for mergers with M tot = 2.4 Mq 
the radius of a nonrotating NS with 1.35 Mq represents 
a better quantity to characterize EoSs in the sense that 
a tighter correlation is found (Fig. [22]) . Correspondingly, 
for M to t = 3.0 Mq the radius of a 1.8 Mq NS is preferable 
as EoS characterizing feature (Fig. |2"4"]) . These trends are 
understandable from the density regimes probed by the 
merger remnant. 

The effects of the initial mass ratio of the binaries are 
only partially addressed in this paper, and a full explo- 
ration is left for the future. Moreover, further studies 
should involve an even larger set of microphysical EoSs 
especially models with a consistent description of ther- 
mal effects. Also, a treatment in full general relativity 
is desirable, preferentially with the inclusion of magnetic 
fields and neutrino effects, to consolidate the scalings pre- 
sented in this work. Finally, a detailed examination of 
the possibility to measure the peak frequency is needed to 
establish the practical use of the described concept based 
on the discovered correlations. For instance, it should be 
explored to which extent the simultaneous observation 
of electromagnetic emission (a gamma-ray burst and/or 
optical transient) and the knowledge of the preceding in- 
spiral signal can support the GW analysis. Moreover, the 
prospects of detector networks should be investigated. 
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Appendix A: Comparison to fully relativistic 
calculations 

In this appendix we explore the reliability of using the 
conformal flatness approximation (see Sect. [TTJ) instead 
of solving the full Einstein equations. To this end we 
compare the GW emission characterized by the ampli- 
tude and the frequency extracted from our calculations 
to results which are available in the literature focusing on 
the postmerger phase as the main subject of this work. 
From earlier works it is known that the quadrupole for- 
mula underestimates the GW amplitude by about 30- 
40 %, whereas the frequency is determined in agreement 
with more sophisticated GW extraction methods [35| . In 
order to judge on the quality of the conformal flatness 
approach we compare the peak frequencies obtained for 
the APR, SLy4, FPS, H3, H4 and Shen EoSs to the val- 
ues found in fully relativistic simulations [H, HI, HI, H3 ■ 
Note that in the calculations of [H[ an ideal-gas index 
r t h = 1.357 was used, for which reason we perform addi- 
tional runs with this value for a better comparison. The 
values of / poa k are listed in Tab. IVIII together with the 
results of fully relativistic studies. If the peak frequencies 
are not explicitly given in the mentioned references, we 
extract them from the provided spectra. 

Note that it is important to take into account differ- 
ences in the exact implementation of the EoSs: In [l8l.l35j 
analytical fits to the microphysical models have been 
used, while piecewise polytropes have been employed 
in HI]. For the FPS, Sly4 and APR models we utilize 
tables provided for these EoSs. The corresponding dif- 
ferences in the implementation of the EoSs lead to differ- 
ences in the stellar structure and thus affect also the peak 
frequency. For instance, for a NS with 1.6 Mq we find 
that the fit to the Sly4 EoS yields a stellar radius which 
is 160 m smaller than the one obtained with our EoS ta- 
ble. Therefore, in Tab. IVIII we recalibrate our values for 
/peak according to the differences found in the TOV so- 
lutions for 1.6 Mq NSs for the various implementations 
of the EoSs. A difference in the radius is translated into 
a difference in the frequency by using the fit formula de- 
rived in Sect. IV B] When such a correction is applied, the 
original value is also given in parentheses. One observes 
a very good quantitative agreement with maximum de- 
viations in /peak of at most a few per cent. Note that 
one cannot identify a systematic over- or underestima- 
tion of the peak frequency by our simulations based on 
the conformal flatness approximation. 

The comparison of the H3 EoS deserves a further com- 
ment. In the fully relativistic simulation of [H| the DRO 
collapses about 5 ms after the merging, whereas in our 
calculation the collapse occurs only after 12 ms. In line 
with the reasoning in Sects. Ivl and IVIll that an early col- 
lapse broadens the peak structure, one observes a double- 
peak in the spectrum of [H| while we obtain a somewhat 
sharper peak (FWHM=255 Hz). The peak at the lower 
frequency in |38| coincides very well with / pea k found 
in our simulation, which is why we conclude that the 



broadening of the peak is connected to the onset of the 
collapse during the GW emission phase. We stress that 
the determination of the exact collapse timescale is sensi- 
tive to the numerical treatment (and e.g., corresponding 
numerical viscosity leading to angular momentum redis- 
tribution) [l9|, HH ■ Therefore differences in the particular 
case of the H3 EoS are not surprising. (Note that the H3 
EoS is excluded because of its low M max . The latter is 
likely to be responsible for the early collapse.) A similar 
reasoning applies to the 1.35-1.35 Mq merger with the 
Sly4 Eos, for which in (3f| the gravitational collapse oc- 
curs within the first 10 ms after merging whereas in our 
simulation the DRO remains stable until the end of the 
simulation (14 ms after merging). 

Notice also that the uncertainties in / pca k due to differ- 
ent resolutions and initial data can be up to a few 10 Hz 
(see [l9| and Sect. IIVB[) . Furthermore, we note that in 
the calculations of [871 ] an approximate treatment of neu- 
trino effects has been taken into account. These may 
affect the exact value of / pea k, but the detailed effect of 
neutrinos on / pca k is not clear from the results of [87^ . 

Our approach predicts a prompt collapse in close agree- 
ment with fully relativistic simulations, e.g. for the bi- 
nary with two 1.5 Mq NSs described by the APR EoS 
(see Tab. IVII|) . In [35( a prompt collapse for the mod- 
els FPS1313 and Slyl414 has been reported, whereas we 
observe for these models the formation of a DRO. How- 
ever, a slight increase in the total binary mass (M to t = 
2.64 M Q for FPS, M tot = 2.9 M for Sly4) is already suf- 
ficient to lead to a prompt collapse. Note that the tran- 
sition between prompt and delayed collapse is crucially 
determined by the EoS. Thus, the observed discrepancies 
in the collapse behavior may be explained by the slight 
differences in the implementations of the EoSs, where 
larger deviations are found for FPS and Sly4 in compar- 
ison to the APR EoS (compare original and recalibrated 
values of / pca k in Tab. IVIl[ ). 

Despite the differences and the remaining uncertain- 
ties it is reassuring that codes with different degrees 
of sophistication yield quantitatively very similar re- 
sults. Note that not only the gravity treatment (con- 
formal flatness approximation versus full general relativ- 
ity) differs between our implementation and the codes 
used in [lH, |35|, [H, H3] , but also different hydrodynamics 
schemes are employed (SPH versus grid-based methods), 
and GWs are extracted by different methods. Without 
more detailed exploration it is therefore unclear to which 
extent on the one hand the different treatment of grav- 
ity and on the other hand the different hydrodynamics 
schemes, the different numerical implementations and the 
different resolutions cause the observed differences in the 
peak frequency. After all, the good match between fully 
relativistic calculations and the conformal flatness ap- 
proximation is not really astonishing in the light of the 
findings of [88| , where agreement of both methods in de- 
termining the frequency of the fundamental quadrupole 
mode has been reported in the case of isolated NSs. 
Moreover, the conformal flatness approximation is known 
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TABLE VII: Models for a comparison between our work and 
fully relativistic calculations. The different setups are char- 
acterized by the employed EoS and the binary components 
Mi and M2 . The second column provides the ideal-gas index 
used for an approximate treatment of thermal effects, "full" 
refers to a consistent description of temperature effects. In 
the third column peak frequencies obtained in our simula- 
tions are listed where we correct for slight differences in the 
EoS implementation (see text). Values in parentheses are the 
uncorrected peak frequencies. /peak aturc denotes the peak fre- 
quencies found in the literature with the corresponding refer- 
ence given in the last column. Models marked by "X" lead to 
a delayed collapse within the first 10 milliseconds after merg- 
ing, "prompt BH" indicates the prompt formation of a BH 
in the given model. 



EoS Mi - M 2 


r th 


/peak [kHz] 


/pcS atm ' c [kHz] 


reference 


APR 1.35-1.35 


2 


3.45 (3.46) 


3.35 


[18] 


APR 1.35-1.35 


1.357 


3.69 (3.68) 


3.58 


[38] 


H3 1.35-1.35 


1.357 


2.67 


2.69-3.00 X 


[38] 


H4 1.35-1.35 


1.357 


2.56 


2.57 


[38] 


APR 1.5-1.5 


2 


prompt BH 


prompt BH 


[18] 


Sly4 1.2-1.2 


2 


3.01 (2.94) 


3.1 


[35] 


Sly4 1.35-1.35 


2 


3.39 (3.32) 


3.6 X 


[35] 


FPS 1.2-1.2 


2 


3.47 (3.41) 


3.5 X 


[35] 


Shen 1.35-1.35 


full 


2.19 


2.18 


[87] 


Shen 1.5-1.5 


full 


2.39 


2.30 


[87] 


Shen 1.6-1.6 


full 


2.57 


2.49 X 


f87] 



to yield exact results in spherical symmetry. 

For assessing the reliability of the GW amplitudes ex- 
tracted in our calculations we directly compare the am- 
plitudes for 1.35-1.35 Mq mergers with the APR, the 
H3, and the H4 EoSs with r th = 1-357 to the results 
shown in [38|, where the GW signal has been computed 
by means of the complex Weyl scalar. This compari- 
son confirms earlier findings that GW amplitudes are un- 
derestimated by the quadrupole formula by about 30-40 
percent [35| • In the ringdown phase we find even slightly 
larger differences. 



Appendix B: Effect of the thermal ideal-gas 
component 

As mentioned in Sects. llTl and Mil not all EoSs include 
temperature effects consistently, for which reason these 
models are supplemented with an ideal-gas component 
to mimic thermal effects. This description requires the 
specification of the corresponding ideal- gas index T t h- A 
choice in the range from 1.5 to 2 seems reasonable fl9| . 
Throughout this paper T t h is fixed to 2, which can be 
considered as suitable approximation but may be slightly 
too high for EoSs yielding compact NSs fl9j . 

To illustrate the uncertainties connected with the ap- 
proximate treatment of thermal effects, Fig. 12"?! shows the 



TABLE VIII: Simulations for different EoSs with varied ideal- 
gas index P t h for mimicking thermal effects, "full" in the 
second column provides results with a consistent temperature 
treatment. All models refer to a binary setup with two stars 
of 1.35 Mq. EoSs providing the full temperature dependence 
(class i) are restricted to zero temperature and neutrinoless 
beta-equilibrium. / poa k denotes the dominant GW frequency 
in the postmerger phase. Note that models which lead to a 
prompt collapse with r t h = 2 are not simulated and not listed 
here, "prompt BH" indicates the prompt formation of a BH 
in the given model. 



EoS 


r t .h 


/peak [kHz] 


GS1 


full 


2.10 


GS1 


2 


2.09 


GS1 


1.5 


2.18 


LS375 


full 


2.40 


LS375 


2 


2.34 


LS375 


1.5 


2.41 


Shen 


full 


2.19 


Shen 


2 


2.16 


Shen 


1.5 


2.26 


GS2 


full 


2.53 


GS2 


2 


2.52 


GS2 


1.5 


2.65 


LS220 


full 


2.89 


LS220 


2 


2.73 


LS220 


1.5 


2.92 


LS180 


full 


3.26 


LS180 


2 


3.08 


LS180 


1.5 


3.34 


eosL 


2 


1.84 


eosL 


1.5 


1.95 


eosO 


2 


2.66 


eosO 


1.5 


2.77 


eosUU 


2 


3.50 


eosUU 


1.5 


3.70 


SKA 


2 


2.64 


SKA 


1.5 


2.81 


APR 


2 


3.46 


APR 


1.5 


3.63 


BurgioNN 


2 


3.20 


BurgioNN 


1.5 


3.33 


Sly4 


2 


3.32 


Sly4 


1.5 


3.53 


Glendnh3 


2 


2.33 


Glendnh3 


1.5 


2.51 


BBB2 


2 


3.73 


BBB2 


1.5 


prompt BH 


eosC 


2 


3.33 


eosC 


1.5 


3.83 
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FIG. 27: Peak frequency of the postmerger GW emission ver- 
sus the radius of a nonrotating NS with 1.6 Mq for fully 
microphysical EoSs (classes (i) and (ii)). Red symbols show 
the results for simulations including thermal effects consis- 
tently. Squares mark the frequencies for computations with 
an approximate treatment of thermal effects using an ideal- 
gas index of Tth = 2. Diamonds correspond to runs with 
r t h = 1.5. Small symbols indicate EoSs which are excluded 
by the observation of the 1.97 Mq pulsar. 



results of additional simulations of 1.35-1.35 Mq merg- 
ers with T th = 1.5 (sec Tab. IVIIIj) . The plot displays 
the peak frequency as a function of the radius of a non- 
rotating NS with 1.6 Mq for all purely microphysical 
EoSs. For all temperature dependent EoSs of our sample 
(Shen, LS180, LS220, LS375, GS1 and GS2) we include 
additional calculations, where the EoSs are restricted to 
the zero-temperature slice and are supplemented with the 
ideal-gas ansatz for thermal effects. For these EoSs we 
construct barotropic relations by assuming in addition 
neutrinoless beta-equilibrium and run simulations with 
r t h = 2 and T t h = 1-5. Peak frequencies obtained from 
these simulations using T t h = 2 are shown with squares, 
while diamonds denote the results of r t h = 1-5 runs. 
Red crosses indicate the peak frequency when the full 
temperature dependence of the EoS is considered. Small 
symbols correspond to EoSs which are excluded by the 
1.97 Mq pulsar. For these models based on class- (i) EoSs 
the r t h = 1-5 results appear at higher frequency com- 
pared to the fully consistent model, and the T t h = 2 
data points are located at lower frequencies. An ideal- 
gas index T t h = 2 appears to be a good choice except for 
the LS EoSs, for which a slightly lower r t h seems more 
suitable. Note that the LS180 EoS is excluded. 

In addition, we simulate for all EoSs of class (ii) (fully 
microphysical EoSs at zero temperature) the symmet- 
ric merger with M to t = 2.7 Mq using T t h = 1.5, which 
should be considered a firm lower bound for T t h [19| . The 
corresponding peak frequencies (diamonds in Fig. I27[) oc- 
cur at higher values, but the differences to the r t h = 2 
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FIG. 28: Peak frequency of the postmerger GW emission ver- 
sus the radius of a nonrotating NS with 1.6 Mq for fully mi- 
crophysical EoSs (classes (i) and (ii)). Symbols have the same 
meaning as in Fig. [9] The data points for EoSs of class (ii) 
(black symbols) are an average of the results for calculations 
with Fth = 1.5 and F t h = 2. 



runs (squares) seem to increase slightly with smaller Ri.e- 
Note that for the BBB2 EoS a prompt collapse to a BH 
occurs for T t h = 1-5. The data points at Ri.e = 13.98 km 
belong to the Glendnh3 EoS whose behavior one should 
take with a grain of salt because of its qualitatively dif- 
ferent mass-radius relation as argued in Sect.fV] For the 
accepted models the largest difference between peak fre- 
quencies extracted from r t h = 2 and r t h = 1-5 simula- 
tions is 210 Hz for the Sly4 EoS. However, the T t h = 1.5 
data points should be considered as safe upper limits for 
the peak frequencies. Because of the lack of better al- 
ternatives we favor the use of Tth = 2 in this work until 
more EoSs including a consistent description of temper- 
ature effects become available. 

In Fig. [28] we display the relation between / pea k and 
Ri.e for the case that the peak frequencies for class (ii) 
EoSs (black symbols) are averaged between the runs with 
r t h = 1.5 and T t h = 2 approximately estimating the out- 
come of a choice of T t h = 1.75. Also in this case a very 
tight correlation is found, but note that in comparison 
to Fig. [Tl] the EoSs excluded by the pulsar mass mea- 
surement of [lj| scatter more strongly from the accepted 
models. Ignoring the result with the Glcndnh3 EoS one 
finds that the relation for the accepted models obeys ap- 
proximately a power law with a maximum deviation from 
this functional relation of 126 meters. The scatter may 
be reduced by choosing a slightly different radius of non- 
rotating NSs for characterizing the EoSs. 

Moreover, we note that in comparison to nucleonic 
EoSs the exact choice of the thermal ideal-gas index T t h 
for strange quark matter is less relevant. Neglecting thcr- 
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mal effects for example in a merger of two stars with 
1.35 Mq described by the MIT60 EoS, decreases the peak 
frequency by only 80 Hz (37j . 
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